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PREFACE 


Ludwig  Prandtl,  concerned  about  the  degree  of  validity  of  certain 
assumptions  in  his  lifting  line  theory,  proposed  for  W.  Kinner  (1937)  the 
problem  of  the  circular  wing  in  it a  exact  formulation.  Kinner  determined 
overall  lift  and  moment  but  did  not  solve  the  span  loading  problem;  his 
results  are  infinite  series  which  appear  to  diverge  at  the  wing  tips. 
Unfortunately,  the  answer  had  to  lie  in  the  tip  solution,  since  the  key 
problem  is  the  transition  between  the  lifting  line  concept  (reap,  its  modern 
extension,  the  standard  collocation  analysis)  and  the  Max  Munk  concept  of 
the  slender  wing. 

The  gap  remained  until  Jordan  (1971)  showed  that  the  apparent  diver¬ 
gence  can  be  overcome.  It  turned  out  that,  contrary  to  both  the  lifting 
line  result  and  the  slender  wing  result,  the  correct  span  loading  has  a 
logarithmic  component  (which,  incidentally,  causes  a  vigorous  impulse  to  the 
vortex  trail  roll-up). 

In  the  present  work,  the  1971  analysis  is  clarified,  extended,  and 
brought  to  its  logical  conclusion.  The  complete  asymptotic  description 
of  the  pressure  singularity  at  the  wing  tip  is  constructed. 

The  wing  tip  pressure  singularity  is  composed  of  elements  of  progres¬ 
sively  increasing  order.  Its  leading  part  is  general  in  the  sense  that  the 
amplitude  ratios  of  its  components  are  fixed  numbers.  Different  solutions 
are  different  only  in  the  specific  part  which  consists  of  components  of  higher 
order.  Our  solution  is  complete  in  the  sense  that  the  leading  general  part 
is  completely  described,  and  the  orders  of  the  leading  components  of  the 
specific  part  also  are  given.  A  practical  consequence  of  these  analytical 
results  is  the  ease  with  which  one  may  calculate  specific  numerical  solutions 
routinely  and  to  very  high  accuracy;  only  a  few  numbers  are  required  to  de¬ 
scribe  such  a  solution  completely. 

Although  the  main  analysis  of  this  paper  deals  with  the  circular  wing 
in  incompressible  flow,  the  form  of  the  general  solution  for  wings  of  arbi¬ 
trary  aspect  ratio  is  deducted  readily  from  the  results  and  allows  ,  for 
example,  discussion  of  the  transition  between  the  finite  aspect  ratio  wing 
and  the  6 Under  wing  model. 


ili 


LIST  OF  CONTENTS 


Page 


ABSTRACT 

PREFACE 

i 

iii 

I. 

INTRODUCTION 

I 

II. 

SURVEY  OF  THE  ANALYTICAL  PROBLEM  FOR  THE  GIF  ’ '  LAR  WING 

2 

III. 

REFERENCE  FORMULAS 

6 

IV.; 

THE  OPERATOR  SQ 

10 

V. 

THE  DOWNWASH  CONDITIONS 

15 

VI. 

THE  MISSING  SET  E£ 

18 

VII. 

INTERPRETATION 

22 

VIII. 

THE  OUTER  PROBLEM 

24 

IX. 

THE  PLANAR  SOLUTION 

26 

X. 

ADDITIONAL  SOLUTIONS 

29 

XI. 

GENERALIZATION  OF  THE  SOLUTION 

30 

XII. 

CONCLUSIONS 

32 

REFERENCES 

4  Tables  and  8  Figures 

33 

APPENDIX  A  ELEMENTARY  SPAN  LOADING  FUNCTIONS  (4  Figures) 

A1 

APPENDIX  B  SUMMATION  FORMULAS 

B1 

APPENDIX  C  ELEMENTARY  SUMS  Cn 

Ujh 

Cl 

APPENDIX  D  PRESSURE  DISTRIBUTIONS  (4  Figures) 

D1 

iv 


L 


1 


I.  INTRODUCTION 

This  paper  is  concerned  with  the  theory  of  thin  lifting  surfaces  of 
finite  span  ("wings")  in  linear  subsonic  flow.  The  distribution  w(£,y)  of 
the  downwash  over  the  wing  surface  is  given.  The  given  downwash  implies  a 
field  of  flow  disturbances  and,  consequently,  a  distribution  p(t,y)  of 
the  aerodynamic  pressure  force  which  acts  on  the  wing.  The  specific  problem 
at  hand  is  to  find  P  <J  ,y).  Because  of  analytical  difficulties  which  arise 
in  particular  at  the  wing  tips,  no  complete  solution  for  any  problem  of 
this  type  had  been  available. 

A  specific  solution,  namely,  the  solution  for  the  planar  circular 
wing  in  incompressib le  flow,  has  been  presented  recently!!) .  The  circular 
wing  represents  a  truly  three-dimensional  (3-D)  problem  in  that  it  is  well 
removed  from  either  of  the  two  limit  cases  for  which  solutions  are  known: 
the  two-dimensional  (2-D)  case  of  a  uniform  wing  having  infinite  aspect 
ratio  and  the  slender  wing  limit  (zero  aspect  ratio).  Also,  the  circular 
wing  represents  a  case  of  a  wing  having  parabolic  wing  tips.  The  parabolic 
wing  tip  is  of  particular  interest  because  it  fits  the  usual  assumption 
(which  stems  from  lifting  line  theory  as  well  as  from  Mink's  minimum  drag 
lemna)  that  the  span  loading  over  any  finite  wing  would  be  of  an  elliptic 
type  (i.e.,  would  fit  an  ellipse  at  the  wing  tips).  From  a  technical 
point  of  view,  the  fit  should  minimize  the  tip  drag;  also,  because  of  the 
fit,  one  would  expect  the  analytical  problem  to  be  less  difficult  than  for 
any  other  wing  tip  shape. 

Nevertheless  there  is  a  reason  why  the  analytical  problem  of  the 
parabolic  wing  tip  has  remained  an  intriguing  challenge  for  many  years; 
the  assumed  fit  does  not;  in  fact  exist.  The  actual  span  loading  is  not  of 
the  elliptic  type.  A  corrective  (logarithmic)  term  was  derived  in  Ref,  (1).®* 

In  the  present  paper,  the  analysis  of  Ref,  (1)  is  generalized  and 
is  brought  to  a  logical  conclusion.  Already  the  corrective  term  of  Ref.  (1) 
is  a  general  term  in  the  sense  that  it  arises  not  from  the  specific  downwash 
w  =  const,  of  the  planar  wing  but  from  a  common  property  of  all  technically 
meaningful  downwash  distributions  w,  A  second  such  general  term  is  derived 
in  the  present  paper.  Not  only  do  such  general  terms  describe  the  struc¬ 
ture  of  the  general  solution,  but  to  know  them  reduces  the  numerical  work 
required  to  calculate  specific  solutions.  (It  also  greatly  reduces  the 
number  of  numerical  values  which  one  has  to  list  to  describe  a  given  spe¬ 
cific  solution.)  On  the  other  hand,  no  such  practical  benefits  arise  from 


(J  The  corrective  term  is  in  conflict  with  a  standard  assumption  of  collo¬ 
cation  analyses.  This  discrepancy  is  not  of  major  concern  where  the  goal 
is  to  find  the  pressure  load  over  the  main  part  of  the  wing,  but  it  markedly 
affects  the  roll-up  mechanism  of  the  vortex  trail. (2) 


knowing  general  terms  which  are  of  higher  order  than  the  differences 
between  specific  solutions.  It  will  be  shown  that  with  the  determination 
of  the  second  corrective  general  term  this  practical  limit  iB  reached. 

A  sequence  of  specific  solutions  has  been  calculated.  These  provide 
insights  of  technical  interest  into  the  mechanism  of  lifting  surfaces. 

Although  the  present  analysis  deals  specifically  with  the  circular 
wing  in  incompressible  flow,  the  conclusions  regarding  the  structure  of  the 
general  solution  can  be  formulated  to  apply  to  arbitrary  wing  planforms 
with  parabolic  wing  tips  in  subsonic  flow.  This  is  briefly  discussed. 

Mathematically,  the  analysis  deals  with  a  type  of  infinite  series 
(a  combination  of  power  series  and  Fourier  series)  which  does  not  seem  to 
have  been  explored  extensively  before.  A  by-product  of  our  investigation 
are  two  sets  of  formulas  of  general  mathematical  interest.  These  formulas 
connect  the  sums  of  infinite  progressions  involving  the -function  (and 
thus  the  logarithm)  to  the  Riemann  ^'function.  They  are  derived  in 
Appendix  B  and  seem  to  be  the  first  formulas  of  this  kind. 


II.  SURVEY  OF  THE  ANALYTICAL  PROBLEM  FOR  THE  CIRCULAR  WING 

In  order  to  be  able  to  give  a  survey  of  the  analytical  problem  at 
hand,  we  start  this  section  by  describing  briefly  one  aspect  of  the  results 
of  Ref.  (1).  In  the  subsequent  section,  some  additional  formulas  of  Ref. 

(1)  are  listed  (sometimes  in  modified  form  resp.  supplemented)  for  easy 
reference.  For  brevity,  we  refer  to  Eq.  (n)  of  Ref,  (1)  as  (l.n).  For 
further  details  and  for  proofs  see  Ref.  (1), 

We  deal  with  the  circular  wing  in  incompressible  flow.  It  was 
assumed  in  Ref.  (1)  that  the  given  downwash  is  symmetric  and  does  not  con¬ 
tain  wing  camber.  Since  no  essential  aspect  of  the  analytical  problem 
is  lost  by  these  restrictions,  we  retain  them  here.  Thus,  using  the  wing 
coordinates  Fig.  1,  we  have  (1.4) 

w(£,y)  5  w(y)  =  w(-y)  (i) 

The  pressure  distribution  p(Jf,y)  which  arises  from  w(y)  can  be  described 
(see  Eqs.  (9)  and  (14)  below)  by  an  infinite  set  of  coefficients 
H  =  0,1, 2, 3...  .  These  coefficients  are  Legendre  coefficients  of  the 
span  load  distribution,  see  Eq.  (9a)  below.  The  leading  coefficient,  Cq, 
can  be  determined  from  the  subset  C^u,  H  =  1,2,3..,  by  means  of  the 
condition  (1.15) 


oo 

C2K.  =  0 


K=0 


(2) 


3 


The  subset  remains  to  be  determined.  We  call  it  the  "solving  set"  for  the 
given  downwash  w(y) . 

Ordinarily,  the  coefficient  Cq  is  positive,  since  it  measures  the 
total  lift:  Cq  3  CjVB,  scc  (l*6a).  Hence  the  elements  of  the  solving  set 
are  ordinarily  negative.  They  can  be  written  in  the  form 


-  C2K^~2  +  RH 


<H*  1)  (3) 


Here  a2  is  a  constant.  The  remainder  has  to  converge  faster  than  the 
leading  term: 


R»  =  o(>C2) 


as 


K' 


*<x> 


(3a), 


-  2 

Thus  the  C2^  converge  as  X  . 

We  describe  the  pressure  p(£,y)  by  the  npn-dimensional  pressure 
function  p  -  (1-Jf2)  ap/q  where  the  factor  eliminates  the  trivial 

singularities  of  half-orders  at  Che  leading  edge  (l.e.)  and  the  trailing 
edge  (t.e.).  The  function  p  is  shown  for  the  planar  case  (vi(y)sl)  in 
Fig.  2.  This  relief  diagram  exhibits  an  interesting  type  of  tip  singu¬ 
larity,  a  singularity  which  is  characteristic  for  the  pressure  distribution 
near  any  parabolic  wing  tip.  This  singularity  is  produced  by  the  leading 


term  of  order  in  Eq.  (3);  its  amplitude  p-^  (1)  is  related  to 

constant  a2  and  to  the  tip  value  of  the  local  lift  coefficient  C^(y)  by 


the 


6  (1)  -  87Ta„  =  2Cj(  1) 

l.e.  2  * 


(4) 


see  (1.34). 

The  relations  Eq.  (4)  are  valid  for  any  arbitrary  set  of  numbers 
C2J.  which  obey  Eqs.  (3,3a),  resp.  for  the  pressure  function  p  which  one 
would  calculate  from  such  an  arbitrary  set.  This  arbitrary  pressure 
function  p  is  finite  everywhere.  There  is  thus  a  temptation  to  assume 
that  any  such  set  C^y  would  represent  the  solving  set  for  a  technically 
meaningful  lifting  surface  problem.*  However,  this  is  by  no  means  the  case. 


*  T^e  pressure  distribution  p  which  belongs  to  p  has  a  singularity  of  order 
(1+1  *")^  along  the  l.e.,  but  this  singularity  is  accepted  as  "technically 
meaningful"  in  linear  theory. 


4 


One  recognizes  this  as  one  calculates  the  induced  downwash  w^(y)  of  p  over 
the  wing  surface.  One  finds  that  w^(y)  will  in  general  diverge  toward  the 
wing  tip: 


w^y)-*  t  oo  as  y-»  1-0 

Such  divergence  is  not  acceptable:  the  w^(y)  of  the  solving  set  is  supposed 
to  equal  the  given  downwash  w(y),  and  any  technically  meaningful  given 
downwash  w(y)  is  finite  everywhere. 

In  Ref.  (1)  the  fact  that  w^(y)  has  to  be  finite  was  used  to  deter¬ 
mine  the  leading  term  of  the  remainder  R*.  Incorporating  this  term,  and 
making  use  of  Eq.  (4),  we  re-write  Eq,  (3)  as 


(  1) 


(5) 


The  new  remainder  R^  here  differs,  of  course,  from  that  of  Eq,  (3);  now 
o(hT^).  Referring  to  Eq.  (5),  we  can  now  describe  the  structure  of 


R>_ 

the  solving  set  as  follows:  it  consists  of  a  known  general  term  (the 

term  in  [  J-brackets  divided  by  4TT)  with  an  unknown  amplitude  C^(l),  and 
a  remainder  R„. 


Between  Eq.  (3)  and  Eq.  (5),  the  infinite  set  H.  ^  has  been  trans¬ 
ferred  from  the  remainder  into  the  general  term.  The  first  question  to 
which  the  present  paper  is  addressed  is  this:  i3  it  possible  to  continue 
this  process,  transferring  step-by-step  further  infinite  sets  (hopefully 
of  successively  increasing  order  of  convergence)? 


There  are  indeed  further  general  conditions  which  are  fulfilled,  in 
practice,  by  all  given  downwash  distributions  w(y) .  Significant  for  our 
purpose  is  the  fact  that  all  the  derivatives  of  w(y)  are  finite  (or  zero) 
at  the  tip.  For  each  successive  derivative,  this  statement  translates 
into  an  analytical  condition  of  increasing  severity  for  &2VL  anc*  should 
yield  one  (or  several)  additional  general  sets  of  successively  increasing 
order.  This  expectation  is  in  agreement  with  the  fact  that  such  sets 
describe  asymptotic  properties  of  the  tip  singularity  in  p  which  are  of 
increasing  order  of  convergence  to  zero  as  the  tip  is  approached  (see  Ref.  1, 
Table  I). 


One  is  interested  to  learn  more  details  of  the  general  term  because 
this  should  reduce  the  numerical  work  involved  in  calculating,  for  a 
given  downwash  w(y),  the  remainder  set  and  the  amplitude  C^(l).  On  the 


other  hand,  one  has  to  expect  to  reach  a  point  beyond  which  to  proceed 
would  have  little  practical  purpose.  Namely,  R^  must  contain,  in  addition 
to  all  the  Identifiable  general  sets,  a  specific  signature  which  describes 
the  given  downwash.  Fran  the  expectation  that  this  specific  signature  will 
describe  the  overall  span  load  distribution  rather  than  its  asymptotic  tip 
behavior,  one  concludes  that  it  will  be  expressed  more  by  the  leading 
elements  (H  -  1,2,3...)  of  the  set  rather  than  by  its  tail  behavior  (the 
behavior  of  R^asH"*®).  Nevertheless,  the  specific  signature  has  to  con¬ 
verge  to  some  order,  say  to  the  order  0S,  as  H-*oo.  There  would  be  little 
advantage  in  knowing  details  of  the  general  term  which  converge  more  rapidly 
than  0_. 

S 


The  second  problem  which  has  to  be  addressed  thus  is  to  determine 
the  order  0g.  The  specific  signature  we  define  as  follows:  we  select  as 
reference  downwash  the  planar  downwash  wQ(y)  a  1.  For  a  given  downwash 
w(y)  4  wQ(y),  we  define  its  "void"  downwash  wv(y)  by 


wv(y)  =  w(y)  -  XwQ(y) 


(6) 


with  the  factor  X  defined  by  the  condition  that  the  tip  value  v(l)  of 
the  "void"  span  loading  becomes  zero: 


c,,v(1>  -  <U(1>  -  “c,o<«  -  0 
The  solving  set  for  wv(y)  is  then,  due  to  Eq,  (5), 

"  C2X,  v  “  RX  "  XRX,o  s  RX,v 


(6a) 


(7) 


This  solving  set,  which  is  by  definition  void  of  all  identifiable  general 
terms,  we  use  as  the  specific  signature  of  w(y). 

We  anticipate  here  briefly  the  answers  to  both  questions.  In  (1.18) 
it  had  been  assumed  tentatively  that  the  remainder  of  Eq.  (5)  could  be 
written  in  the  "rational"  form 


oo 

R*  =  ^  (ar/*,r) 

r»4 


(8) 


However,  later  numerical  evaluation  of  the  planar  solution,  Fig.  2,  indicated 
that  Eq.  (8)  was  incomplete:  step-by-step  determination  of  the  constants  ar 
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did  not  converge  in  a  satisfactory  manner.  Indeed,  a  review  of  the  analysis 
confirmed  that  at  least  one  "missing  set"  has  to  be  added  to  Eq.  (8).  In 
the  analysis  of  the  present  paper  this  missing  set  is  determined.  It  is 
O(tt_4log2»0  and  is  the  third  general  set  (it  follows  the  two  general  sets 
which  are  given  in  Eq.  (5)).  To  proceed  further  with  this  analysis  would 
entail  increasing  difficulties.  Fortunately,  the  limit  of  practical 
interest  has  already  been  reached:  the  order  of  the  fourth  general  set 
(appears  to  be  K.-^log^K  and)  is  higher  than  the  order  0g  of  void  sets  which, 
from  numerical  results,  appears  to  be  R"^logjc. 

The  span  loading  distributions  which  belong  to  the  elementary  sets 
mentioned  in  the  preceding  discussion  are  illustrated  in  Appendix  A. 

A  third  aspect  of  the  problem,  also  discussed  in  this  paper,  are 
the  numerical  methods  required  to  obtain  a  specific  solution.  For  the 
reference  ilownwash  wn(y),  one  determines  0/(1)  by  truncating  the  given 
infinite  linear  system  to  N  equations,  matching  the  solution  with  the  known 
general  terms,  and  extrapolating  to  (l/N)  =  0.  For  additional  solutions, 
one  determines  X  (rather  than  C^(l))  in  a  simpler  process.  After  the  re¬ 
mainder  has  been  determined,  its  leading  terms  are  extracted;  one  is 
then  left  with  a  final  remainder  set  which  converges  very  rapidly.  The 
numerical  procedure  is  self-checking  since  the  leading  coefficient  Cq  is 
calculated  directly  and  the  result  must  fulfill  Eq.  (2). 

In  the  present  paper,  three  additional  solutions  are  presented 
numerically.  Thus  one  has  available  a  three-fold  variety  of  solutions  by 
linear  superposition. 


III.  REFERENCE  FORMULAS 


Th 

is  (1.5) 


^  n  L 

e  formula  for  the  non-dimensional  pressure  function  p  =  (l-Jz)'ip/q 


p  =  - 


4  r  -22.  ™  oo  +1  7 

,,  2y-i  C2Mr"  cos  2*0  +  c2A+lr  Sin(2\+1)0  j 

U-y  [g,0  A=0  J 


(9) 


The  local  lift  coefficient  is  (1.6) 


c.(y)  -  .  i  r 

0--y2)  *  **o 


oo 

y  C2*P2R(y) 


(9a) 
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and  the  position  of  the  local  center  of  pressure  is  (1.7) 


i*(y>  "  "  c^yT  <2A+1>C2A+lP2A(y) 


The  P2ll(y)  and  Pq^(y)  are  systems  of  orthogonal  polynomials.  The  former 
are  the  Legendre  polynomials,  the  latter  their  derivatives  (Gegenbauer  or 
ultraspherical  polynomials),  normalized  such  that 


P2>t(±1>  =  *  1 


(all  K,A) 


P2X^  ~  (X+1)(2A+1)  dyp2^+1 


The  two  sets  of  unknown  coefficients  in  Eq.  (9),  the  C21t  and  the 
CjJ^+i,  are  fully  interdependent.  We  introduce  the  abbreviation 

A  a  A  +  h 


Then  (1.13)  reads 


J0  TT 


I  ^  C2A-H  . 


A=0  * 


.2  ^^2i±l 

*  A-0  K2-  A2 


We  have  available  also  the  reversed  form  (1.19)  of  Eq.  (13) 


_2_  ^  d? 
C2A+l  *  7TX  &  * 2 


H=1  *i  -  A 


(Note  that  the  system  Eqs.  (13)  fulfills  Eq.  (2)  for  any  arbitrary  set 

c2A+l-> 
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Usually,  the  span  loading  C^(y)  is  of  more  immediate  concern  than 
the  center  of  pressure  Jf^  (y)  or  the  complete  pressure  function  p.  Hence 
the  set  C2H  would  appear  to  be  of  more  iumedlate  interest  than  the  set 
^2A+1'  However,  Eq.  (9a)  takes  the  form  (0/0)  as  y-^1;  the  important 
tip  value  Cf(l)  is  thus  not  immediately  given  by  the  set  C^. 

For  this  tip  value  we  have  the  formula  (1.22b) 


oo  m 

Ml)  -  8  SrAC2A+1  (15) 

A  =0 


This  formula  has  two  interesting  consequences.  One,  consider  an  arbitrary 
set  RH  which  obeys  Eq.  (3a)  (for  instance,  any  finite  set  R*,  is  eligible). 
According  to  Eq.  (4),  this  set  Redoes  not  contribute  to  (>4(1)  (i.e.,  the 
tip  value  of  its  pressure  function  p  is  zero).  From  R4,  calculate  the  set 
which  corresponds  to  R^  by  means  of  Eq.  (14).  From  Eq.  (15)  follows 


00 

2^  Rk  -  0  (15a) 

A  =0  A 


for  all  such  sets  R^. 

Two,  using  Eq.  (15)  in  Eq.  (9b)  to  determine  the  tip  value  of  the 
pressure  center  function,  we  obtain,  noting  Eq.  (10) 

^X<1>  -  -  T/4  (16) 

for  any  pressure  distribution  (at  least  so  if  C^(l)  f  0).  This  result  is 
in  agreement  with  the  prediction  of  slender  wing  theory  for  the  planar 
elliptic  wing.  As  an  illustration,  the  pressure  center  function  Jf^(y)  for 
the  planar  circular  wing  is  shown  in  Fig.  3.  (The  planar  span  loading  is 
shown  in  Fig.  A3). 

The  given  downwash  w(y)  enters  the  analysis  in  the  form  of  its  set 
of  Gegenbauer  coefficients  ws  (1.39)® 

00 

w(y)  =  2Z  (s+D(2s+l)wsP2g(y)  (17) 

s=0 


0  We  use  here  the  notation  w,  in  place  of  the  notation  w  of  Ref,  l. 

s 
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For  example,  for  the  planar  wing  (i.e.,  for  our  reference  downwash  wQ(y)) 
we  have 


w0  -  1  ;  s  0  if  w(y)  s  wQ(y)  s  1  (17a) 

The  aet  of  given  dovmwash  coefficients  ws  and  the  unknown  set  C2K  are  inter¬ 
related  by  the  Infinite  linear  system 


w 


s 


1 

+ - 

K+s+i,- 


(18) 


see  (1.12)  with  (1.37).  The  notation 


s  s  s  + 

is  here  used.  It  corresponds  to  the  notation  X»  Eq.  (12);  in  fact,  the  two 
indices,  s  and  A ,  play  rather  similar  roles  in  the  analysis  (as  we  will  see 
shortly) . 

The  system  Eq.  (18)  is  a  convenient  basis  for  our  analysis  because 
of  the  relative  simplicity  of  its  matrix.  It  is  directly  equivalent  to  the 
(less  convenient)  original  Kinner  system,  Eq.  (60)  of  Ref.  3.  The  alternate 
system 


w 


s 


cAC2A+1 


(19) 


see  (1.14),  relates  the  given  set  wg  to  the  unknown  set  This  system 

is  preferable  for  numerical  work.  It  converges  well  because  its  matrix 
differs  little  from  the  unit  matrix  (see  Table  Illb  of  Ref,  1);  on  the 
other  hand,  the  are  too  complicated  for  convenient  analytical  treatment 
(see  Table  II  of  Ref.  1). 
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IV.  THE  OPERATOR  Sr 


Using  operator  notation,  we  can  write  Eq.  (18)  as 


ws*S<<W;  C2R*"S  <wa>  (20) 

Our  problem  would  be  solved  if  we  could  find  an  explicit  form  for  the 
reversed  operator  S-^.  Realistically,  we  cannot  expect  to  find  this  explicit 
form  in  spite  of  the  relatively  simple  form  of  S.  However,  we  can  determine 
the  salient  features  of  S"*. 

The  operator  S  would  be  still  simpler  if  the  sumnand  %  in  the  denomi¬ 
nator  of  the  second  term  of  its  matrix  would  be  zero.  Indeed,  when  <K+s) 
is  large,  this  4  would  appear  to  be  negligible.  Actually,  as  we  will  sec, 
all  the  complications  which  make  our  analytical  problem  so  Interesting  are 
caused  by  this 

Separating  out  the  complication  we  have  for  Eq.  (18)  the  form 


-  ^2  'C2K  1  • 
Ws  =  2s  27  "  2'  Z2  +  1  2 
K»0  A  -s4  *  M  i 


oo  r  1  i  -j 

=0  *•  H+5+.V  K+s  J 


s0<C2h)  +  S1<C2*> 

Elimination  of  Cq  by  means  of  Eq.  (2)  yields 


oo  -u  C 

S0<C2H>  “  <2/S)  ^  -j-rp 

11=1 


(21a) 


,  °o  r  i 

sl<c2*>  - -  - 

2  k=l  L s(k+s) 


(k+s)  (i+^>)  (K+S+’i) 


]  "  C2K 


(21b) 


Of  these  two  operators,  Sq  is  dominant  in  Eq.  (21);  is  of  higher  order. 
With  Sq  we  are  already  familiar,  and  we  also  know  already  its  reversed 
operator  Sq  :  we  can  write  Eq.  (14)  as 

C2V+1  s  (1/n“)S0(C2M)  ;  C2K  55  ^*s0  (c2^+l^  (22) 

and  can  read  Sq1  directly  by  comparing  Eq.  (22)  with  Eq.  (13). 
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It  follows  that  we  possess  already  a  first  order  approximate 
solution,  namely,  the  solution  which  arises  when  we  set  5  0.  The  set 
C2^  can  then  be  calculated  from  the  given  set  wR  by  means  of  Eq.  (20), 
and  the  set  is  equivalent  with  the  set  wR  since  both  sets  arise  from 

^2H  ^  means  of  the  same  operator  (setting  a  0  is  equivalent  with 
replacing  the  matrix  cjj  in  Eq.  (19)  by  the  unit  matrix  from  which,  in  fact, 
cjjj  does  not  differ  much).  The  brief  table  which  follows  illustrates  the 
degree  of  approximation  which  is  reached  in  this  manner  for  our  reference 
case,  Eq.  (17a).  The  first  line  of  the  table  gives  the  formula  for  each 
coefficient,  the  second  and  third  lines  give  the  (rounded-off)  numerical 
values,  approximative  and  correct: 


coefficient 

c0 

C1 

-c2 

c3 

-C4 

C5 

-C6 

approx.  / 

2/H2 

0.2026 

i/n 

0.3183 

4/3  ff2 

0.1351 

0 

0 

4/15H2 

0.0270 

0 

0 

4/351 r2 

0.0116 

correct 

0.2238 

0.3496 

0.1446 

0.0041 

0.0307 

0.0014 

0.0135 

(23) 

One  sees  that  the  error  made  by  setting  s  0  is  of  the  order  10%. ® 

Because  of  the  thus  illustrated  importance  of  the  operator  S0,  we 
set  out  to  investigate  its  properties  in  some  detail.  Since  our  results 
will  be  accurate  relations  between  the  sets  C^+l  and  C2K  but  only  approxi¬ 
mate  relations  between  ws  and  C2^,  we  will  write  our  analysis  in  terms  of 
^2A+1  (however^ we  will  later  use  these  results  in  discussing  the  relations 
between  ws  and  C2tt) . 

We  use  the  notation  E*  for  any  one  of  the  elementary  sets  of  wh  . 
the  solving  set  -Co*  may  be  composed.  We  do  not  stipulate  that  has  to 
be  a  rational  set  (i.c.,  E^  =  H”r)  but  assume,  on  the  basis  of  Eq.  (3) 
and  of  numerical  experience,  that  EK  is  a  progression  in  H.  which  converges 
to  zero  smoothly  as  K  -►  00 .  We  use  the  integer  r  to  describe  the  range  of 
convergence  between  x-r  (included)  and  (excluded),  writing 

EH=  K'rg|t  (r  &  2)  (24) 


with  gR  =  o(*0  but  not  o(l) .  If  Ea  is  a  rational  set  we  have  5  1;  we  may 
also  have,  for  example,  gK  =  logn*  or  g*  =  k  2. 

The  contribution  to  which  arises  from  by  means  of  Eq.  (14) 

we  d  mote  by  E’  c^us 


0  See  also  Appendix  D, 
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IT  c 


2M-1,E 


oo 

-  (2/K)  X 


*-l  <H2-I2)^r‘2 


(25) 


The  sum  on  Che  right  we  take  apart  by  means  of  the  identity 


(2fl) 

<K2-*2)Hr"2 


-  (l/A1^1)  S  fl  +  (-)r'‘'J  (|)l> 
*  -2  L 


+  <1/Ar) 


(26) 


Terms  with  i>  arise  if  r  i  4.  Summation  of  these  terms  in  Eq.  (25)  does  not 
pose  any  difficulty®  and  lead3  to  rational  elementary  sets  JT~(t-#+l) 

®2A+1J  with  (r-iH-1)  odd  and  1  3.  More  interesting  Is  the  final  term  on 

the  right  of  Eq.  (26).  For  its  contribution  we  use  the  following  abbreviation: 


T  C2A+1,E  "  •••  +  (x^r)  °j,,X 


Thus 


<8«> 


2  /-i-  -  i 
ti  U-JT  h 


(27) 


From  the  point  of  view  of  analytical  difficulties,  the  operator  0 r  (gH) 
forms  the  nucleus  of  the  operator  Sq.  ’ 

The  details  of  the  reversed  operator  Sq^-  are  similar  to  those  of  Sq. 
Corresponding  to  Eq.  (24),  we  write  Ex  for  the  elementary  sets  of  anc* 

E,  -X-tgA  (tfe3)  (28) 


<2  If  g„»  1,  summation  leads  to  the  Riemann  J -function.  If  gg  is  more 
complicated,  it  may  be  necessary  to  use  Eq.  (B2). 
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From  Eqs.  (13)  and  (22) 


"  T C2lt;E 


2^ _ 

A=0  (K2-A^)Jt':L 


(**)) 


(29) 


An  identity  similar  to  Eq,  (26)  is  used  to  transform  this  sum.  Its  terms 
with  P  produce  ra;ional  elementary  sets  in  C2x.  There  remains  a  final 
operator  which  is  rimilar  to  the  operator  0^  ^(g4) : 

-  Tc2i^e  °  —  +  (i/xF)0^^) 


with 


The  two  O-operators,  Eqs.  (27)  and  (30),  are  readily  executed  in 
the  case  of  rational  elementary  sets  (i.e.,  gk  =  1  resp.  gx  5  1).  They 
then  lead  into  the  realm  of  the  psi  (digamma)  function.  One  findB 


with 


and 


(  4L*t  )  f  even  > 

=  1  o  I lf  ‘  “  f ^  i 

(31) 

L, -1  +  i  +  i+..+  2>ii.i^(,+i)  - Y (%>j 

^(logXf^)  +  log  2  +  1/48*2  +  .... 

Ola) 

fCl/A))  r  even  •» 

V*(I)  *  Uj  I if  r  is  ?  odd } 

(32) 

=  +  ( ^/^A)  »  log  2 

(32a) 

where 
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The  significant  observation  here  is  that  the  0-operators  can  produce 
logarithmic  elements  from  rational  elements  Eg  reap.  Ej^. 

If  we  apply  the  preceding  results  to  the  tentative  assumption,  Eq.  (8), 
that  only  rational  elements  H."r  occur  in  62^  we  find  that  the  corresponding 
elementary  contributions  ^ave  t^'e  following  form: 

**TC2X+1,E  "  <A3^3>  +  <Af/A5)  +  (A^/A7)  +  ...  (33) 

The  factors  At  are  given  by  the  following  matrix: 


A3 

A5 

a£ 

r»2 

0 

0 

0 

3 

•41* 

0 

0 

0 

4 

•2J(2) 

1* 

0 

0 

5 

-2*(3) 

-4l£ 

0 

0 

6 

-2y(4) 

-2J(2) 

0 

7 

-2f(5) 

-210) 

-4£ 

0 

As  expected,  this  matrix  contains  logarithmic  elements.  It  follows  a  simple 
law  and  is  easily  extended;  also,  it  is  easily  reversed. 

•  -1 

Conversely,  no  logari chmic  elements  can  arise  if  the  operator  Sq 

is  applied  to  Eq.  (33)  to  regain  the  elements  of  C2g.  This  explains,  see 
Eq.  (31),  why  only  terms  with  t  odd  occur  in  Eq.  (33J.  Interesting  is 
the  observation  that  the  operator  Oj  K(g^)  even  eliminates  the  logarithm 
from  the  elements  with  Ljy  .  To  this  observation  we  will  return  later. 

The  main  result  of  the  preceding  investigation  of  the  operators 
So  and  So  is  that  their  nuclei  are  the  two  0-operators.  A  stimulus  for 
investigating  the  latter  further  will  arise  when  we  consider  the  complete 
relations  between  the  sets  C2g  and  ws .  Before  we  turn  to  this  task,  let 
us  utilize  Eq.  (33a)  to  take  a  look  at  the  actual  solutions,  the  sets 
C2*  and  C2A+1j  ^or  our  reference  case  Eq.  (17a). 

Let  Ug  be  the  negative  of  the  general  term  in  the  n  -bracket  of 
Eq.  (5): 

G0  -  +  £f<2)  -  *<3>/4j 
Gg  -  -  [*-2  -  H'3/4  ] 


(«*0) 


(34a) 
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In  Fig.  4,  the  ratio  C^/G*  is  plotted  over  M-  For  K-+  ao,  this  ratio 
should  converge  towards  C,<l)/4tT*  0. 12677.  Indeed,  after  an  initial 
disturbance  around  M«*l,  convergence  is  both  rapid  and  smooth.. 

From  Eq.  (33),  the  set  G^  which  arises  from  G^  is 

Ga  -  [l  +  L*  J/l TA3  (all  A )  (34b) 

The  ratio  C^+l^A  is  also  plotted  in  Fig.  4.  There  is  again  smooth  con¬ 
vergence,  but  convergence  is  new  slow.  Indeed,  all  the  elements  X,~r  con¬ 
tribute  elements  0(A”^)  to  Cjji+i,  see  Eq.  (33a).  Thus,  while  the  components 
A"3logA  should  have  corresponding  amplitudes  in  C^-f^  and  GA,  already  the 
components  0(A"^)  are  likely  to  be  different.  Slowness  of  convergence  is 
an  obvious  consequence. 

A  final  remark:  the  formulas  of  this  section  can  be  used  to  derive 
certain  sumnation  formulas  of  mathematical  interest.  This  is  discussed 
in  Append  ijc  B. 


V.  THE  DOWHWASH  CONDITIONS 

In  view  of  the  normalization  Eq.  (10),  it  follows  from  Eq.  (17)  that 
the  downwash  coefficients  w8  must  converge  o(s‘^)  if  the  tip  value  w(l)  is 
to  be  finite.  If  one  differentiates  Eq.  (17)  once,  one  finds  that  the 
derivatives  of  the  P2S  are  of  order  s^  at  the  tip;  hence  in  order  that  w'(l) 
be  finite,  the  wB  must  converge  o(s_5).  Further  differentiations  impose 
further  convergence  conditions  of  increasing  severity.  However,  as  already 
indicated,  there  will  be  no  need  for  us  to  drive  our  analysis  beyond 
0(s"^). 


We  have  to  use  Eq.  (20)  in  order  to  translate  the  convergence  con¬ 
ditions  for  ws  Into  conditions  for  C2 ^  At  first,  we  make  again  the 
tentative  assumption  that  C2*  is  composed  of  rational  sets  E^,  Eqs.  (6), 
(8).  For  these  sets,  we  know  already  approximate  elements  wg  g,  namely, 
the  elements  7I‘C2A+1J  E  of  Eqs.  (33),  (33a).  (We  have  to  replace  A  by  s, 
of  course.)  To  these  we  have  to  add  as  corrections  the  contributions  of 
the  operator  Sj.  As  we  do  this,  the  regular  form  of  Eq.  (33a)  is  lost, 
and  there  is  then  no  longer  an  advantage  in  using  the  abbreviations  s  and 
L*  .  Therefore  we  write  our  results  in  terms  of  s  and  log  s: 


w 


s,E 


_  1 

2f, 


00 


j\ log  s  +  +  C*  J/b> 


(s  *  0) 


(35) 


16 


It  turns  out  that  the  matrix  of  coefficients  Bj  is  (about)  a  tri¬ 
angular  matrix,  with  zeros  below  a  zigzag  line  along  its  diagonal.  The 
matrix  Cj  is  a  full  matrix.  Of  the  two,  we  will  use  explicitly  only  the 
following  coefficients: 


aI 

<S 

K 

r=2 

-1 

+2.5 

2.25 

r=3 

-4 

-J<2) 

7.5 

ri4 

0 

-|>t<r-2>  +3(r-!)J 

0 

(35a) 


Both  matrices,  Bj  and  Cj,  must  have  their  basis  in  the  matrix  A*>, 

Eq.  (33a).  Specifically  the  mechanism  which  produces  Bj  is  as  follows: 
this  matrix  must  arise  from  the  log- terms  -4Lt  in  Eq.  (33a).  For 
example,  the  coefficient  B^  corresponds  directly  to  Aj.  From  this,  contri¬ 
butions  to  ...  arise  in  two  ways.  One,  from  developing  the  factor, 

X'3  in  Eq.  (33),  8"^  in  the  case  of  wg  E,  in  terms  of  powers  s"^, ...  . 
Two,  from  the  operator  S^. 


The  mechanism  is  the  same  if  r=2,  in  principle  at  least,  even  though 
there  is  no  log-term  for  r=2  in  Eq.  (33a).  The  point  is  that  in  principle 
there  is  a  log-term  A^.  Its  factor  is  zero  because  now  r  is  even  in  Eq.  (32). 
But  Eq.  (32)  does  not  apply  to  the  operator  Sp  In  consequence,  creates 
non-zero  higher  order  coefficients  ...  .  x 

The  mechanism  for  all  r  odd  is  similar  to  that  for  r=3,  and  for  r 
even  similar  to  that  for  r=2.  The  result  is  an  almost  triangular  matrix  Br. 
The  first  non-zero  coefficients  on  successive  lines  r  are  B^,b3,B^,B|,B^ 


By  a  similar  mechanism,  a  full  matrix  arises  from  the  constant 
coefficients  in  Eq.  (33a).  Because  of  this,  it  is  readily  possible  to 
fulfill  successively  the  conditions  that  the  rational  terms  s"3,s“^,  ... 
in  ws  should  have  the  factor  zero,  Combine  Eqs.  (5),  (8)  to  read 

"St*  ZI<ar/x^  50  that  ws  =  21  (arwS)E>  (36) 

r=2  r=2 


d  This  is  as  described  in  the  text  which  precedes  (1.42).  Regrettably, 
the  manner  in  which  the  order  of  the  remainder  terms  are  written  in  both 
(1.42)  and  (2.21)  is  incorrect.  Both  equations  should  be  corrected  to 
agree  with  Eq.  (35a). 
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One  simply  hag  to  use  a  sufficient  number  of  elementary  sets  iCv  and  has 
to  determine  their  amplitudes  such  that  the  conditions 


arC^  =  0  (36a) 

are  fulfilled  for  *11  required  values  & . 

Conversely,  it  is  not  possible  to  extract  more  detailed  conditions 
for  the  ar  from  Eq.  (36a).  Hence,  in  what  follows  we  can  disregard  the 
rational  terms  as  neither  causing  complications  nor  providing  useful  informa¬ 
tion. 

Specific  conditions  arise  from  the  log-terms,  however.  First,  the 
factor  to  the  terms  e'^log  s  is  zero®,  from  Eq.  (35a),  if  and  only  if 

a3  =  "a2/4  (37) 

This  result  is  already  incorporated  into  Eq.  (5). 

In  order  to  proceed  beyond  this  result  of  Ref.  1,  we  have  to  eliminate 
next  the  term  s'^log  s.  Here  we  encounter  an  obstacle.  Because  of  Eq.  (37) 
we  have  now,  due  to  Eq.  (35a), 


ws  =  (3/16)  a2  s'^log  s  +  ...  (a  f  0)  (38) 

This  term  is  different  from  zero  because,  in  general,  a2  =  C^(l) /4JV  0.  No 
further  set  is  available  in  Eq.  (35a)  to  cancel  this  term. 

The  inescapable  conclusion  is  that  Eq.  (38)  requires  an  additional 
elementary  set  in  C^,  a  set  that  is  "foreign"  to  the  rational  sets  K"r 
which  are  represented  in  Eq.  (35a).  Applied  to  this  "missing"  set,  the 
operator  S  must  produce  a  contribution  which  cancels  the  leading  term  in 
Eq.  (38)  but  must  not  contain  lower  order  terms  (disregarding  lower  order 
rational  terms) . 


0  We  omit  the^*  of  Eq.  (35)  for  simplicity.  If  we  wish,  we  may  consider 
as  part  of  the  rational  term. 
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VI.  THE  HISSING  SET  E* 

ft# 

The  task  of  finding  the  missing  elementary  set  would  appear  to 
imply  that  one  has  to  try  various  likely  sets,  apply  the  operator  5,  and 
hope  to  obtain  a  result  suitable  to  cancel  the  term  Eq.  (38) .  Since  for 
transcendental  sets  the  required  summation  formulas  are  not  usually 
available,  the  task  as  thus  described  would  appear  to  be  formidable. 
Fortunately,  we  have  available  an  argument  that  reduces  the  same  task  to 
a  direct  form  which,  furthermore,  requires  only  the  performance  of  analyti¬ 
cal  integrations. 

In  view  of  Eqs.  (26/7),  we  write  the  missing  set  as 


E*  -  l'* 


s* 


(39) 


The  operator  Sq  produces  the  following  contribution  to 


irC2A+l,E*  "  -27g<2>/*3  +  VJW/A4 


(39a) 


where 


OO  gy 


Assume  that  we  can  find  gK  such  that  the  result  of  the  0-operator  in  Eq. 
(39)  produces  log  A  as  leading  term.®  Then  the  corresponding  contribution 
of  the  operator  will  be  of  the  general  order  A’3;  this  contribution 
is  thus  not  here  of  interest.  Hence  we  can  use  the  same  O-operator  for  ws; 
the  leading  terms  here  are  thus  found  to  be 


w 


s.E* 


•  -2(jg(2)  +  %  5gO)J  /s3  +  s-4log  s  +  ... 


(39b) 


0  The  logarithmic  sets  and  Lj  which  w^  used  earlier  were  defined  for 
all K  and  all  A  j  in  particular,  Lq-0  and  Lq  =  \  -  log  2.  Just  as  we  used 
log  s  instead  already  in  Eq.  (35),  we  use  henceforth  also  log  A  even  though 
log  0  does  not  exist.  This  formal  procedure  simplifies  the  analysis,  and 
it  is  justified  because  we  are  concerned  only  with  the  asymptotic  behavior 
as  s  resp.M  resp.A  oo  .  It  is  readily  shown  that  in  this  respect  the 
formal  analysis  leads  to  the  correct  result. 
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The  log-term  in  Eq.  (39b)  can  be  used  to  cancel  the  term  Eq .  (38). 

Further,  since  we  had  to  use  only  the  0-operator  of  Sq,  not  of  S^, 
we  can  reverse  the  process.  Again  disregarding  rational  sets,  we  have  for 
the  missing  set  from  Eq.  (30),  since  gB  *  gt/7 T, 


8*  =  °s,k.(1o2  s>/7?2  (t  =  4) 


(40) 


Equation  (40)  puts  the  task  of  finding  the  missing  set  into  a  direct 
form.  This  form  still  involves  inconvenient  infinite  summations  over  pro¬ 
gressions  which  involve  log  s.  However,  again  because  we  are  only  inter¬ 
ested  in  the  leading  term,  we  can  deal  with  these  summations  by  the  method 
of  corresponding  integrals.  (Essentially,  this  means  using  only  the  lead¬ 
ing  term  on  the  right  of  Eq.  (B2).) 

With  t=4,  we  have  from  Eqs.  (30),  (40)  formally 


irV  -  2  (~ 

s=0  'a 


— i—  -  — !—  )  log  S 
s-K.  s+*J 


(41) 


This  sum,  after  some  rearrangement,  including  folding  of  its  middle  term, 
becomes 


fcl  ,  M-l  , 

25:  -  21  -1- 


log 


K+S 


OO 


=0  s+^  s-0  a+^  «"8-l  s  =  n  S+% 


-  Z  los  (l-(«/s)2) 


The  corresponding  integral  of  the  second  of  these  three  sums  is 


1  2 

I  log  dx=  /  log  1+u  du  .  JL  +  00T1)  (42a) 

if  «-%-x  x  (2JI-1)"1  1_u  u  4 


The  integral  which  corresponds  to  the  last  sum  is 


f  log(l- (K/x)  )  =  /  log(lV)  _ du _ 

x+%  0  u(l+u/2n) 


JL-  +  oOc"1) 

12 


(42b) 
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Neither  contribution  is  here  of  interest.  The  leading  contribution  arises 
from  the  first  sum: 


*-l  i-h 

2  J  ISJLJZ  dx  ■*  2  (  log(u-^)  =  log^jt  +  0(1)  (42c) 

0  ^  u 


As  we  had  to  expect,  the  leading  contribution  is  not  a  rational  set.  Also, 
note  that  this  contribution  arises  as  s-*ie;  clearly  our  formal  procedure 
was  justified. 

We  can  now  list  the  sets  (leading  terms  only)  which  belong  to  E*  in 
ws  (from  Eq.  (39b)),  C2n  and  C^+i: 


ws,E* 

“C2X,E* 

C2A+1,E* 

l  o-K-  a 

lo&i* 

l  og  A 

s* 

i r2*4 

*A4 

(43) 


This  list  completes  the  task  which  was  set  by  Eq.  (39).  We  use  wg  E*  to 
cancel  the  undesirable  term  on  the  right  of  Eq.  (38)  and  arrive  at' a  more 
detailed  form  of  Eq.  (5) : 


16VZHT  J  4  W 


(*>1) 


(44) 


To  the  new  set  E*  in  Eq.  (44)  belongs,  because  of  Eq.  (43),  a  set 
in  C^+ij  which  is  not  contained  in  Eq.  (33a).  This  term  was  indeed  found 
in  a  numerical  analysis  of  the  set  Vfa+l’  Fig.  4. 

Equation  (44)  essentially  completes  the  task  of  describing  the  struc¬ 
ture  of  the  solving  set  C2j^.  It  disproves  the  tentative  assumption  Eq.  (8) 
and  implies  that  to  proceed  beyond  Eq.  (44)  would  require  a  more  involved 
analysis.  In  particular,  in  using  the  method  of  corresponding  integrals 
for  E*  we  neglected  the  higher  order  terms  which  we  would  than  need. 
Fortunately,  as  already  mentioned,  there  is  no  practical  need  to  drive  the 
analysis  any  further. 

We  round  off  this  discussion  of  the  sets  E*  iry  showing  that  the 
relation  between  them,  Eq.  (43),  can  also  be  derived  by  applying  the  method 
of  corresponding  integrals  in  reversed  order. 
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I 


From  Eqs.  (2.7)  and  (43) 

IT^s  E*  +  <1/b4>  S  ( -i-  -  log2*  (45) 

*  K"l'X-5  X+5  J 

By  rearrangement,  the  sum  becomes 

lgK2(2§)  +  JL  f  log2(s+*0  -  log2(s-K+l)l 
i  lt-1  t- 


+  2SEI  —  r  log2(H+s) 
*.=8+2  L 


log2(K-e-l)J 


The  single  term  is  of  no  interest  and  can  be  dropped.  The  corresponding 
integrals  for  the  sums  are 


]  £log2(s-bc)  -  log2(s-x)  J  ^  +  J  £log2(x+s) 


log2(x-s)7  ^ 
X 


Setting  x=us  in  the  first,  ux=s  in  the.  second  and  combining,  we  arrive  at 
the  following  integral 


J^2  log  5  +  log (~  -  u)J 


log 


1+u  du 


after  the  limits  of  the  two  integrals  have  been  replaced  by  their  values 
for  s-^oo.  We  have  here  the  integral  Eq.  (42a),  multiplied  by  4  log  s, 
and  a  further  integral  which  is  0(1)  and  not  here  of  interest.  Inserting 
in  Eq.  (45),  we  confirm  Eq.  (43). 
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VII.  INTERPRETATION 

Since  the  solving  sets  C2*,  converge  with  increasing  slowness  as  r, 
increases,  a  lengthy  tabulation  would  be  required  should  one  wish  to  com' 
municate  a  given  numerical  solution  without  knowing  the  structure  of  C2R, 
Eq.  (44).  Because  of  Eq.  (44),  one  has  to  list  only  the  number  0^(1) 
plus  the  much  faster  converging  set  R^.  We  refer  to  the  problem  of  deter¬ 
mining  the  structure  as  the  "inner  problem"  and  to  the  problem  of  calculat¬ 
ing  C^(l)  as  the  "outer  problem".  With  the  latter  we  will  deal  in  the 
next  section  and  conclude  the  discussion  of  the  former  in  the  present 
section. 

First  we  collate  our  results,  Eqs.  (31),  (33a)  reversed,  (43)  for 
the  operator  Og  For  this  purpose  let 

ws  =  s^gg  and  -C2ll  =  K^g*/?2 


Again  only  the  leading  terms  of  each  set  are  listed: 


Es 

®  M,  t  even 

Sr,  t  odd 

1 

2  log  R 

0 

log  s 

log2* 

-1/2 

(46) 


In  general,  the  are  0(lognR)  and  are  thus  within  the  convergence  range 
defined  for  g-sets.  In  all  these  cases,  neither  O-operator  alters  the 
general  order  of  convergence,  i.e.,  r“t.  The  two  exceptions  are  borderline 
cases,  gs=l  with  t  odd,  and  g*=l  with  r  even,  see  Eq.  (32). 

The  operator  Og  ^  increases  the  power  of  the  log-function  by  one  if 
t  is  even  and  decreases  it  by  one  if  t  is  odd.  Its  reversed  operator 
0^,5  increases  the  power  when  r  is  odd  and  decreases  it  when  r  is  even. 

This  observation  we  extrapolate  to  speculate,  in  a  heuristic  manner, 
about  further  st  uctural  details  beyond  those  given  in  Eq.  (44).  Sets 
K'^logH  and  r,"^  produce  only  rational  sets  s~*  (if  we  disregard  t  1  5), 
There  is  thus  no  obvious  reason  to  expect  these  sets  in  the  general  term 
of  Rather,  if  we  accept  the  scheme  Eq.  (35a)  up  to  0(s"^),  then 

(again  because  the  matrix  is  a  full  matrix)  these  two  sets  are  the  first 
sets  for  which  there  is  no  detailed  general  condition  for  their  amplitudes. 
To  this  observation  corresponds  the  (numerical)  result  that  these  two  sets 
are  the  leading  sets  in  the  remainder  R^  of  Eq.  (44). 
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It  follows  that  the  next  general  term  should  be  of  the  general  order 
r*5.  Note  that  E*  produces,  besides  s"^log  s,  also  a  potential  term 
(having  the  factor  zero)  s-4log3s.  The  latter  spawns  a  term  s"^log^s  via 
the  operator  Sj_;  to  cancel  this  term,  a  term  R"^log^  is  required.  This 
then  would  appear  to  be  the  fourth  general  term. 

Whether  or  not  this  speculative  result  is  correct,  the  point  is  clear 
that  complications  would  increase  rapidly.  With  these,  we  are  not  going  to 
involve  ourselves  any  further.  Instead  we  rewrite  Eq.  (44)  somewhat  for 
the  purpose  of  representing  numerical  solutions. 

A  modification  is  that  we  combine  the  first  two  general  terms  in  a 
manner  which  produces  a  simpler  contribution  to  wfl: 


C2X,E 

^E 

1  1 

H2  4H3 

£0  +  1(2)] 

s'3  +  1.5  s“4log  s  ... 

4 

X(4*+l) 

10  s“3 

+  1.5  s-4log  s  . . . 

The  difference  between  the  two  forms  of  C2j,  £  is  0(*"4)  and  is  thus  not  of 
concern  in  the  general  term  (it  does  of  course  affect  R^) . 

From  the  remainder  Rj^we  split  off  its  two  leading  terms.  Thus 


•  C2H.  " 


r, 


H(4H+1) 


3  log2*  I  _1 

i6  r2*4  J  47 r  +  *4 


(47) 


To  represent  a  given  numerical  solution,  we  have  to  list  the  three  constants 
C*(l),  b^  and  plus  the  final  remainder  R^.  We  will  find  that  the  latter 
converges  very  rapidly. 
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VIII.  THE  OUTER  PROBLEM 

Even  though  Eq.  (47)  gives  more  details  of  the  structure  of  the 
general  term  than  Eq.  (5),  it  would  still  be  difficult  to  extract  its 
amplitude  0^(1)  from  a  given  (and  necessarily  truncated)  numerical  set  02^. 
The  use  of  Eq,  (15)  and  thus  of  the  set  appears  preferable. 

We  mentioned  already  above  at  the  end  of  Section  II  that  for  a  given 
wing  one  has  to  determine  C^(l)  only  for  the  reference  downwash  wg(y)  ■  1. 
In  Ref.  1  the  extrapolation  formula  (1.52)  was  used  for  this  purpose  and 
led  to  the  approximate  amplitude 


0^(1)  ft  1.5931  (48) 

We  show  in  this  section  how  this  value  can  be  improved  by  utilizing  our 
newly  derived  knowledge  about  the  structure  of  the  set 

Write  Eq.  (15)  in  the  form 


C*<1) 


S(N)  +  a(N) 


N 

4  j?  (2A+1)C 

0 


2M1 


00 


+  4  ?E  (2A+1)C 
N+l 


2\+l 


For  a  given  (truncated)  numerical  set  CjA+lf  tlie  sum  ^(N)  is  racily  calcu¬ 
lated.  To  explore  the  structure  of  a(Njy  write 


(2A+l)C2A+i 


DglogA  +  D2logA  +  D3  + 

( 2A+1>  2  (2A+1)3 


(49) 


in  agreement  with  our  earlier  results.  Using  Eq.  (B2),  then 

A(N)  -  f(l  +  log  N)Dfl  +  D,]  H  .  +  . .  . 

(  2N+1)  2 

This  shows  that  the  curve  S(N),  plotted  over  u  ■  (1/N),  has  a  vertical 
tangent  at  u»0.  Its  extrapolation  to  u=0  (in  order  to  obtain  Cj(l)  ■  S(oo)) 
would  be  rather  unreliable. 

The  vertical  tangent  is  eliminated  if  one  uses  the  method  of  Ref.  1, 

writing 


CjCl)  =  s*(N)  +  a*(N) 


(50) 
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with 


S*(N)  -  S(N)  +  4H(2H+l)C21ffl 

Then 


a*(N) 


4N+3  „ 

- o  D0 

(2W1)2 


[dog  N  -  ^)D2  +  D3J 


(2N+1)" 


(50a) 


(50b) 


Tlie  curve  a*(N)  has  the  form  0(u),  and  the  curve  S*(N)  therefore  has  the 
form  Cg(L)  -  0(u) .  Extrapolation  has  now  become  much  more  reliable. 

To  improve  the  procedure  of  Ref.  1,  we  make  use  of  the  fact  that  the 
constants  Dq  and  D2  in  Eq.  (49)  are  directly  related  to  the  unknown  Cj(l); 

D0  *  C^D/TT2  ;  D2  -  -30^(1) /47T 2  (51) 

from  Eqs.  (33), (37)  and  (4)  for  Dq,  Eqs.  (43),  (44)  for  D2.  In  other  words, 
we  match  our  extrapolation  procedure  to  the  known  structure  of  the  set 
C2A+11  (The  constants  Dl  and  D3  are  related  to  the  remainder  R^  rather  than 
to  the  general  term  and  are  hence  unknowns ;  note,  however,  that  does  not 
appear  in  the  leading  terms  of  Eq.  (50b).) 

Inserting  first  only  Dq,  we  write 

c.(l)  -  S*(N) /  [l  -  1  +  a  lo«  M-:-b  +  ...  (52) 

K  »-  T2(2N+1)2J  (2N+1)2 

We.  can  plot  the  first  term  on  the  right,  S*/£J,  and  know  that  its  curve, 
extrapolated  to  u=0,  has  to  have  a  horizontal  tangent  at  its  end  point  C^(l). 

The  steps  which  we  discussed  so  far  are  illustrated  in  Fig.  5.  The 
almost  vertical  line  S*  at  the  left  is  the  interpolation  curve  of  Ref.  1. 
There  are  two  curved  marked  S*/CJ;  these  illustrate  the  remarkable  increase 
in  reliability  and  accuracy  that  is  achieved  by  inserting  the  relation 
between  Cjj(l)  and  Dq.  The  lower  of  the  two  curves  is  constructed  using  the 
values  the  results  of  solving  the  reversed  system  Eq.  (19),  truncated 

at  K  =  N.  The  curve  through  the  points  which  have  thus  been  calculated  for 
a  number  of  values  N  must  reach  <u<«  at  u=0,  but  we  do  not  know  its  shape 
near  its  end  point. 

N 

To  obtain  the  final  (upper)  curve,  the  values  were  extrapolated 

to  N  =  oo  for  each  X-  Through  a  remarkable  coincidence,  it  was  possible  to 
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construct  a  curve,  through  the  thus  calculated  points,  using  only  the  last 
written  out  terra  on  the  right  of  Eq,  (52),  without  allowing  for  any  addi¬ 
tional  terras,  and  using  constant  numbers  a  and  b.  This  curve  is  shown  in 
Fig.  5.  It  connects  all  calculated  points  to  within  drawing  accuracy  (no 
points  were  calculated  beyond  N  =  40)  even  though  the  vertical  scale  of 
Fig.  5  is  highly  stretched.  The  first  noticeable  deviation  between  points 
and  curve  occurs  to  the  right  outside  Fig.  5,  at  N  2  10. 

The  end  value  of  the  new  curve  is  C^(l)  ©  1.5930904.  It  confirms 
the  result,  Eq.  (48),  of  the  less  sophisticated  extrapolation  of  Ref.  1. 
However,  it  is  clear  that  the  new  value  cannot  be  quite  accurate.  One,  we 
cannot  expect  that  the  higher  terms  in  Eq.  (52)  are  entirely  negligible; 
two,  the  number  a  in  Eq.  (52)  should  be  related  to  the  end  result  C^(l) 
via  Eq.  (51).  Making  use  of  point  two  and  allowing  for  point  one  we  found 
Cx(l)  a  1.5930884.  The  smallness  of  the  difference  between  the  two  new 
values  Indicates  the  order  of  accuracy  which  one  achieves  with  this  extrapo¬ 
lative  method.  The  last  decimal  of  the  last  value  is  to  be  considered 
uncertain. 

In  the  tabulation  of  the  planar  solution,  Table  I,  we  will  use  the 
amplitude  number 


a2  -  0.1267740  (Table  I)  (53) 

From  it  follows  0^(1)  *  4lT  a2  ■  1.5930890..  .  This  number  a2  is  treated 
in  Table  I  as  an  exact  number.  It  allows  to  construct  to  8  decimals 
even  though  a2  itself  is  given  only  to  7  decimals.  To  this  point  we  will 
return  later. 


IX.  THE  PLANAR  SOLUTION 

N 

In  the  last  section  we  used  the  truncated  set  C2^+^  which  was  obtained 
by  solving  the  truncated  linear  system  Eq.  (19).  The  corresponding 
set  cJL  is  calculated  by  means  of  Eq,  (13).  Again  the  final  solution  C^has 
to  be  obtained  by  extrapolation.  However,  convergence  as  N-»oo  is  con¬ 
siderably  slov/er  in  the  case  of  C2^  than  in  the  case  of  C^+i,  and  the  re¬ 
quired  extrapolative  procedure  deserves  a  brief  description. 


Two  successive  results,  C2J^  and  C^, 
reasons.  One,  corresponding  elements  C 

new  term  C2tj+i  is  added  in  Eq,  (13)  when 
second  cause  has  by  far  the  larger  effect, 
for  N  large  on  the  basis  of  Eq,  (13)0 


differ  because  of  two  independent 


the 


and  C”  differ.  Two, 

1VU  c£  The 

Taking  this  into  account,  we  write 


0  A  factor  \  has  to  be  added  when  M  =  0.  However,  the  extrapolative  pro¬ 
cedure  is  the  same  for  all  M- 
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+  C 


N-l 

2  K 


JL  2NC‘2N-fl 

T  N2  -X.2 


(1  -*) 


(54) 


The  number  S  is  introduced  here  to  cater  for  the  first  effect.  Numerical 
testing  shows  that  6  is  small  compared  with  1. 

To  simplify  further,  we  replaced  in  Eq.  (54)  by  from  Eq. 

(49)  with  the  higher  constants  D^,  D3,  ...  all  set  equal  to  zero.  Using 
this  simplified  form,  we  calculated  6  for  N  *  37,38,39,40.  We  found  that  6 
changed  little  with  K,  and  changed  even  less  with  N. 

If  one  assumes  that  #  is  sufficiently  independent  of  N,  one  can, 
using  Eq.  (B2) ,  sum  for  each  M  the  to  N=oo  and  obtain  a  tentative  final 
solution  C^.  To  check  the  result,  we  repeated  the  procedure  using  N  =  60. 
Since,  within  our  accuracy  requirement,  the  result  was  the  same,  we  accepted 
it  as  the  final  result. 

We  calculated  the  A  to  12  decimals.  About  2  decimals  are  lost  in 
the  sunmation.^  The  further  calculation  was  done  carrying  10  decimals.  The 
final  result,  R||,  is  tabulated  to  8  decimals  in  Table  I. 

In  Table  I,  the  planar  solution,  presented  in  the  form  of  Eq.  (47), 
is  denoted  by  C2jj  q  to  distinguish  it  from  the  other  solutions.  Tables  II  to 
IV,  The  set  C2lt  g  is  tabulated  to  K=  15,  and  the  remainder  (which  is 

C2{|,o  Plus  the  general  part)  is  tabulated  in  the  second  column. 0 

To  split  off  the  leading  terms  of  R*,  see  Eq.  (47),  one  proceeds  as 
follows.  One  chooses  a  tentative  pivot  point  H  =  p.  The  amplitudes  b^  and 
are  determined  from  the  two  conditions 


-  0  ; 


Condition  b)  arises  from  Eq.  (2)  due  to  the  fact  that  all  the  sets  in  the 
general  part  as  well  as  the  two  leading  sets  of  R^  fulfill  Eq.  (2)  individ¬ 
ually. 

Of  course,  the  goal  is  that  Rr  =  0  for  all  Xfc  p.  If_p  was  chosen 
too  small,  then  this  will  not  be  the  case.  The  tail  sum  of  Rjj  will 
still  be  zero,  but  the  values  R a  beyond  X  =  p  will  describe  a  wavy  curve 
around  the  zero  line.  On  the  other  hand,  if  p  was  chosen  unnecessarily 
large,  significant  decimals  will  have  been  lost,  and  the  constants  b^  and 
will  be  ill-drfined. 

In  Table  I,  RJt<  0.5x10  ®  already  when  X=  15.  Thus  knowledge  of  the 
structure  of  C2#,  Eq.  (47),  has  enabled  us  to  represent  the  slowly  converging 
set  C2n  in  a  relatively  very  brief  table  to  very  high  accuracy. 


P  _ 

b)  2,  -  0  (55) 

0 


0  The  leading  constants  Cq  g  of  the  here  required  sets  C2-  g  are  given  in 
Appendix  C.  ’  * 
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Table  1  can  be  used  to  reconstruct  to  about  the  accuracy  of  R^, 
that  is,  to  about  8  decimals.  In  doing  this,  one  has  to  treat  the  constants 
a2»  b4  aiK*  c4  as  exact  numbers,  even  though  these  constants  are  given  to  7 
decimals  only  and  are  only  approximations  of  their  mathematically  defined 
counterparts . 

The  inaccuracy  Just  mentioned,  the  differences  between  the  listed 
constants  and  their  mathematically  defined  values,  interferes,  to  the 
accuracy  of  Table  I,  neither  with  the  validity  of  Eq.  (2)  nor  with  the 
manner  in  which  Redoes  converge.  Namely,  these  inaccuracies  affecu  columns 
2  and  3  only  as  far  as  the  first  few  values  (for  X.  =  0,1,2..)  are  concerned, 
and  the  sum  of  these  values  is  zero,  for  each  column,  closely  enough  to 
cancel  out  the  discrepancies. 

1= 

The  first  few  values  in  Table  I  are  not  claimed  to  represent  their 
mathematically  exact  images  accurately  to  8  decimals,  but  they  are  accurate 
to  8  decimals  if  the  purpose  Is  to  reconstruct  C2^>  (At  the  higher  values 
of  X,  of  course,  the  inaccuracies  in  question  affect  only  decimals  beyond 
those  listed.) 

In  the  process  of  calculating  the  set  C 2 including  Cq,  no  explicit 
use  was  made  of  Eq.  (2).  The  fagt  that  it  was  possible  to  extract  after¬ 
wards  a  smoothly  converging  set  RK  which  satisfies  Eq.  (2)  serves  as  an 
overall  confirmation  of  the  numerical  procedure. 

It  might  seem  contradictory  that,  while  it  is  relatively  easy  to 
determine  the  set  C2*.  to  very  high  accuracy,  it  is  considerably  more  diffi¬ 
cult  to  determine  0^(1) .  The  point,  already  mentioned,  is  that  C2*  defines 
the  span  loading  C^(y)  well  enough  over  the  inner  part  of  the  wing  but  not 
near  the  wing  tip.  An  illustration  of  the  tip  region  is  Fig.  6.  Curves 
<U.<y)  are  shown  which  have  been  calculated  with  Eq.  (9a)  truncated  at 
N  *  20,25,30,35,40.  All  these  curves  turn  to  +  00  as  y  -#»  1,  where  cor¬ 
rectly,  from  Eqs.  (5)  and  (A4),  (see  (1.47a)) 


-  fi  +  ri  (i-y2)^  log  ...1 
L  16  1_„2  -* 


VD 


(56) 


To  obtain  the  curve  Cj(y),  one  has  to  fair  the  curve  Eq.  (56)  into  the  curve 
for  N-40,  say,  at  about  y*0.98.  The  result  is  shown  in  Fig.  A3. 

Note  that  the  key  which  makes  this  matching  process  possible  lies 
in  the  fact  that  certain  analytical  functions,  Eq.  (A2),  have  Legendre 
coefficients  which  behave  asymptotically  like  the  rational  sets  R'r. 
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X.  ADDITIONAL  SOLUTIONS 

The  principle  which  allows  to  calculate  additional  solutions  quickly 
once  a  reference  solution  has  been  determined  was  already  indicated  in  Eqo. 
(6)  and  (7).  We  have  used  this  principle  to  calculate  the  three  solutions 
^21^1  f°r  downwaeh  distributions 


vL(y)  =  P2i(y)  (i  =  1,2,3) 


(57) 


Together  with  the  reference  planar  solution  C2lt  q*  ^hese  solutions  form  a 
set  of  four  independent  solutions  (and  thus  represent,  by  linear  super¬ 
position,  a  three-fold  multitude  of  solutions). 

N 

jjjjhe  sets  ^  were^calculated  for  N  ■»  60.  By  comparing  the  tail  of 
a  set  £  with  that  of  C  2*  q,  a  preliminary  value  X  was  obtained  and,  using 
this  value,  a  preliminary  void  set  C1^  iv  was  calculated.  Since  for  this  set 
the  constant  Dq  in  ^2/41  is  (almost)  zero  by  definition,  extrapolation  to 
N*oo  is  simplified.  The  resulting  set  is  treated  essentially  like  the  set 
Rg.in  the  reference  case,  except  that  now  there  is,  in  addition  to  b^  and  c^, 
a  third  parameter,  namely,  a  small  adjustment  aX  to  X.  This  parameter  is 
necessary  to  make  the  final  remainder  R^converge  towards  zero  without  waviness. 
In  this  numerical  step,  we  again  carried  10  decimals.  No  difficulty  was 
encountered  when  p  =  25  was  chosen  as  the  pivot  point.  The  results,  again 
rounded  to  8  decimals,  are  listed  in  Tables  II  to  IV. 

A  graphical  presentation  of  the  four  solutions  i  and  their  respective 
downwash  distributions  w.  is  Fig.  7.  All  curves  Cj  ±  form  a  downwash  hook 
near  the  tip  with  a  vertical  tangent  at  the  tip  itself;  however,  in  the  pre¬ 
sentation  of  Fig.  7  this  hook  is  too  small  to  be  recognizable  when  i*3  or 
i*4.  The  corresponding  void  curves  iv  are  shown  in  Fig.  8. 

The  relative  scales  of  ^  and  w^  are  chosen  in  Fig.  7  such  that 
Cf  0  and  wq  have  about  the  same  amplitude.  For  i?*0,  the  amplitudes  of  the 
wa(res  in  Cj  t  arc  much  smaller  than  the  amplitudes  of  the  corresponding  waves 
in  w^.  The  amplitude  ratios  correspond  roughly  to  the  ratio  between  half¬ 
wave  length  and  wing  span  (this  is  what  one  would  expect  from  slender  wing 
theory) . 

An  interesting  observation  is  the  following:  one  may  form  the  ratio 
between  two  integrals  over  the  wing  area,  the  lift  integral  and  the  downwash 
momentum  integral.  For  a  planar  wing  this  ratio  is  %CL.  For  the  void  parts 
of  the  solutions  i  =  1, 2  and  3,  this  ratio  is  found  to  differ  from  that  for 
the  planar  wing  by  less  than  it.  As  one  calculates  the  same  ratio  for  the 
complete  solutions,  Fig.  7,  small  differences  between  large  figures  have  to 
be  formed,  and  the  resulting  ratio  varies  somewhat  more,  between  0.895  for 
i*=0  and  0.821  for  i=3.  The  indication  remains  that,  as  far  as  the  total 
lift  of  a  given  wing  is  concerned,  the  details  of  how  the  total  downwash  is 
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distributed  over  the  wing  span  are  of  relatively  little  relevance.  The 
details  of  this  distribution  tend  to  have  negligible  effect  on  the  outer  flow. 

We  have  discussed  in  this  section  only  the  span  loading  functions 
C/j(y).  The  complete  pressure  distributions  p(j*,y)  of  closely  related 
solutions  are  discussed  in  Appendix  D,  These  provide  further  interesting 
insight  into  the  mechanism  of  lifting  surfaces. 


XI.  GENERALIZATION  OF  THE  SOLUTION 

Pressure  singularities  of  orders  (l.e.)  and  +*?  (t.e.)  arise  due 
to  linearization  at  those  wing  edges  which  cross  the  direction  of  the  undis¬ 
turbed  flow.  These  singularities  are  well  understood.  The  wing  tip  problem 
is  the  problem  of  the  transition,  between  the  two  types  of  singularities, 
where  the  wing  edge  becomes  parallel  to  the  flow.  This  problem  is  consider¬ 
ably  more  involved.  We  were  able  to  resolve  it  for  the  specific  case  of  the 
circular  wing  because  here  a  complete  set  of  analytical  relations  of  relative 
simplicity  are  available:  the  elementary  pressure  solutions  in  Eq.  (9);  the 
reversible  relation  between  the  two  sets  of  coefficients,  those  for  the  span 
loading  and  those  for  the  chordwise  moments,  Eqs.  (13)  &  (14);  the  relation 
between  amplitude  coefficients  and  downwash  coefficients,  Eq.  (18). 

Corresponding  analytical  tools  are  not  available  for  wings  of  arbi¬ 
trary  planform,  and  it  will  not  readily  be  possible  to  determine  the  details 
of  the  solutions  for  such  wings  to  the  extent  that  we  did  in  Eq,  (47)  for 
the  circular  wing.  On  the  other  hand,  it  is  in  the  nature  of  the  lifting 
surface  mechanism  that  the  type  of  pressure  singularity  which  arises  at  a 
local  wing  singularity  is  not  affected  by  the  far  field  from  other  parts  of 
the  wing.  For  the  circular  wing,  we  found  that  the  leading  part  of  the 
pressure  singularity  is  independent  even  of  the  slope  of  the  downwash  dis¬ 
tribution  at  the  tip  and  is  determined  solely  by  the  geometric  tip  shape; 
the  remainder  part  which  describes  a  particular  solution  goes  to  zero  and 
does  so  of  higher  order  as  the  tip  is  approached.  The  essence  of  this 
observation  must  be  valid  for  all  wings  having  parabolic  wing  tips,  and  the 
reference  length  for  the  coordinates  which  describe  the  tip  singularity  must 
be  the  wing  tip  radius  Rj. 

Since  to  bring  this  qualitative  statement  into  a  detailed  quantita¬ 
tive  form  would  require  extensive  further  investigations,  we  confine  our 
present  observations  to  a  few  main  points.  In  particular,  we  do  not  concern 
ourselves  with  details  of  the  pressure  distribution  (which  we  did  not  discuss, 
beyond  Fig.  2  and  Appendix  D,  even  for  the  circular  wing).  The  span  loading 
distribution  is  more  readily  discussed  since  here  we  have  available  certain 
analytical  functions  which  conveniently  describe  the  leading  terms  of  the 
limit  behavior  at  the  tip. 
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The  analytical  functions  in  question  are  those  whose  sets  of 
Legendre  coefficients  converge  like  the  rational  sets  see  Eqs.  (A2,2a). 
Hence  we  can  represent  the  first  two  general  sets  of  Eq.  (44)  as  already 
shown  in  Eq,  (56)  where,  on  the  right,  we  should  add  a  remainder  function 
R(y) .  However,  we  do  not  have  a  corresponding  analytical  function  for  the 
set  nr^log2*.,  A  further  difficulty,  discussed  in  Appendix  A,  is  that  already 
a  set  0(>t"4)  might  have  a  finite  tangent  at  the  tip,  Fig.  Al,  or  a  vertical 
tangent,  Fig,  A2,  To  decide  on  this  point  one  has  to  review  the  complete 
set,  not  only  its  rate  of  tail  convergence®.  To  avoid  further  discussion  of 
such  higher  order  details  we  confine  our  further  attention  to  the  two  terms 
already  given  in  Eq.  (56). 

Let  b  be  the  wing  span,  so  that  u  =  (l-y)b/2Rr  is  the  proper  spanwise 
coordinate  for  describing  the  wing  tip  singularity  asymptotically  near  the 
tip  y=+l..  Here  (l-y^)fc  2u  for  the  circular  wing.  For  general  wing  plan- 
forms,  we  can  hence  write  the  local  lift^(y)  as 


'(y>  =  +  •^((l-y2)b/2RT)%log  ]  VX)  +  R(y)  }  c(y) 


(58) 


We  use  the  wing  chord  c(y)  since  the  curve  A  (y)  tends,  to  a  degree,  to 
imitate  c(y)  over  the  inner  part  of  the  wing.  Considering  two  wings  having 
the  same  tip  radius,  the  difference  between  the  two  chords  is  0(uc(y))  near 
the  tip  and  is  hence  of  higher  order. 

The  complete  lift  function  X{y)  has  a  vertical  tangent,  but  it  differs 
from  the  parabolic  type  0(uz)  of  c(y)  by  the  log- term.  The  remainder 
function  R(y)  of  Eq,  (58),  taken  by  itself,  may  be  expected  to  go  to  zero 
about  like  0(u2) . 


The  generalized  form  of  the  solution,  Eq.  (58),  is  valid  for  any 
aspect  ratio  and  thus  in  particular  for  the  slender  wing.  For  the  elliptic 
wing,  slender  wing  theory  predicts  an  elliptic  span  loading,  that  is,  Eq. 

(58)  without  the  log-term  and  without  R(y) .  The  log-term  is  of  the  order 
u(AR)zlog  u  if  AR  the  aspect  ratio.  It  follows,  as  a  sample  application  of 
Eq.  (58),  that  slender  wing  theory  is  valid  strictly  only  in  the  limit  AR-*0. 


0  The  void  span  loading  curves  Fig.  8  represent  void  functions  R(y)  and 
appear  to  have  finite  tangents.  However,  it  might  be  that  the  small  com¬ 
ponent  0(K"4logn)  in  RK  enforces  a  vertical  tangent;  this  detail  would  not 
be  visible  in  the  scale  of  Fig.  8. 
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XII.  CONCLUSIONS 

In  the  earlier  paper^^  we  had  shown  that  the  apparent  divergence  of 
the  Kinner^)  series  is  eliminated  by  the  physically  necessary  condition  that 
the  infinite  set  of  span  loading  coefficients  C2K.  has  a  zero  sum,  and  that 
the  properties  of  the  sum  of  the  Kinner  series  can  be  investigated  by  means 
of  the  sets  of  Legendre  coefficients  of  certain  analytical  functions.  We  had 
further  shown  that,  in  order  to  have  a  finite  downwash  inside  the  tip,  the  span 
loading  function  must  have  a  logarithmic  component. 

This  earlier  analysis  is  pursued  further  in  the  present  paper  and  is 
brought  to  itB  logical  conclusion.  The  duality  between  the  set  C2*,  and  the 
set  C2X+1  of  moment  coefficients  is  worked  out,  and  the  close  relationship 
between  the  anc*  t*ie  downwash  coefficients  ws  is  formulated  and  is 

utilized. 


It  is  shown  that  either  C-set  has  to  contain  logarithmic  elementary 
sets  in  addition  to  rational  sets,  and  the  asymptotic  description  cf  the 
pressure  singularity  at  the  wing  tip  is  constructed.  Its  elementary  sets 
are  identified  up  to  the  order  of  the  "void"  part  of  the  solution,  the 
latter  being  that  part  of  a  specific  solution  which  distinguishes  it  from 
a  reference  solution.  The  void  part  converges  to  zero  at  the  tip,  essen¬ 
tially  like  u^/2  if  u  the  distance  from  the  tip  referred  to  the  tip 
planfomi  radius  %. 


With  the  structure  of  the  solutions  thus  identified,  it  becomes 
relatively  easy  to  calculate  numerical  solutions  to  very  high  accuracy,  in 
particular  once  a  reference  solution  has  been  determined.  Furthermore, 
only  a  short  tabulation  is  required  to  fully  describe  such  a  solution. 

The  main  analysis  of  this  paper  deals  specifically  with  the  circular 
wing  in  incompressible  flow.  The  form  of  the  general  solution  for  arbi¬ 
trary  wings  with  parabolic  wing  tips  is  briefly  discussed.  This  form  shows 
in  particular  that  the  result  of  slender  wing  theory  is  valid  strictly  only 
in  the  limit  of  zero  aspect  ratio. 


Certain  conclusions  of  direct  engineering  interest  can  be  drawn  from 
the  calculated  sample  solutions  for  the  circular  wing.  One,  the  details 
of  how  the  total  given  downwash  is  distributed  over  the  wing  span  affects 
the  total  lift  relatively  little.  Two,  a  spanwise  wavy  incidence  distribu¬ 
tion  is  well  reproduced  by  corresponding  waves  in  the  pressure  distribution 
along  the  wing  leading  edge,  but  over  the  rear  part  of  the  wing  no  visible 
waviness  in  the  pressure  remains  (assuming,  of  course,  that  the  wavelength 
is  relatively  short) . 


A  by-product  of  the  analysis  are  two  sets  of  formulas  of  general 
mathematical  interest.  These  formulas  connect  the  sums  of  infinite  pro¬ 
gressions  involving  the  y -function  (and  thus  the  logarithm)  to  the  Riemann 
^-function  and  seem  to  be  the  first  formulas  of  this  kind. 
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Table  I.  Planar  Solution  C2^  q 
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1620 

+ 

33 

12 

1138 

+ 

17 

13 

823 

4 

9 

14 

610 

+ 

5 

15 

462 

+ 

3 

16 

357 

+ 

1 

17 

280 

+ 

1 

18 

222 

(Kil8>  0 

Go. 

.  o..  + 

+  c 

/.)  +  R„  -  0 

b4  =  0.0024522 
c4  =  0.2270359 


w2jV(y)  -  P4(y)  -  X  -  <21y4  -  14y2  +  l)/8  -  0.2795066 


2*,2v 


Table  III.  Solution  C 


X. 

C2fc,3 

** 

0 

-.04502880 

+.43750001 

1 

+.03210567 

“  • 

39520727 

2 

+.00938151 

m  a 

3174875 

3 

+.0105  6414 

“  • 

1494365 

4 

~. 00^19578 

+. 

381883 

5 

98681 

+  . 

42559 

6 

37073 

+  . 

10118 

7 

17716 

+  • 

3217 

8 

9684 

+  . 

1211 

9 

5786 

+  . 

510 

10 

3686 

233 

11 

2466 

+  « 

113 

12 

1716 

+  . 

57 

13 

1232 

+. 

30 

14 

908 

+  . 

16 

15 

685 

+  . 

9 

16 

526 

+  . 

5 

17 

411 

+  . 

3 

18 

326 

+  . 

1 

19 

262 

+  . 

1 

20 

213 

(*L20)  0 

C2*,3v  +  ^b4LM.+  c4>  +  **  =  0 


b4  =  -0.0156775 
c4  =  0.3787791 


w3,v(y)  =  Vy)  •  X 

=  <429y6  -  495y4  +  I35y2  -  5)/64  -  0.2068387 


Table  IV.  Solution  C2h,3v 
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Fig.  2.  Relief  of  Pressuie  Function  p  over 
Right  Half  of  Planar  Circular  Wing 
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APPENDIX  A 

ELEMENTARY  SPAN  LOADING  FUNCTIONS 


In  the  analysis  of  this  paper  the  solving  set  forms  an  inter¬ 
mediary  between  3ivcn  downwasli  w  and  required  pressure  distribution  p. 

The  critical  analytical  problem  arises  in  the  relation  between  w  and  C2^ 
and  the  elementary  sets  which  are  used  to  build  up  C2H  are  taylored 
for  this  critical  problem,  Comparatively,  the  relation  between  C2*.  and, 
say,  the  span  loading  (y)  is  a  simple  mathematical  formula,  Eq.  (9a), 
Nevertheless,  this  formula  does  not  give  the  reader  an  immediate  illustra¬ 
tion  of  the  contributions  of  the  sets  Eh  to  C^(y) ,  The  purpose  of  this 
Appendix  is  to  pro\Hde  this  illustration. 

If  occurs  in  the  general  term  of  Eq.  (6),  then  its  contribution 
to  C^Cy)  is 

p  -l  oo 

ci<y)  =- — 5^2Xp2H.<y)  (Al) 

2(l-y^) 2  o 


if  the  normalization  0^(1)  =  1  is  made.  We  use  Eq.  (Al)  as  our  definition 
of  elementary  span  loading  functions. 

In  Ref.  1  elementary  sets  F2^  were  considered  whose  span  loading 
functions  are  simple  analytical  functions,  see  Table  I  of  Ref.  1. 
Abbreviate  by  writing  y  for  (1-y2) then,  for  r  A  2, 


;r-2 


2CX<y>E=F^"  {  _r_ 


yr'2log(2/y) 


if  r  is  even 
if  r  is  odd 


(A2) 


The  first  six  of  these  functions  are  plotted  in  Fig.  Al.  In  standard 
collocation  analyses,  only  the  functions  with  r-2,4,6..  are  used.  Of  these, 
r»2  denotes  the  strictly  elliptic  distribution  (a  straight  horizontal  line 
in  Fig.  Al) .  The  curve  r*4  has  a  finite  tangent  at  the  tip;  the.  subsequent 
curves  r=6, 8..  have  zero  tangents.  On  the  other  hand,  the  correct  curve 
C*(y)  has  a  vertical  tangent  at  the  tip.  Clearly,  one  has  to  include  the 
function  r=3  if  one  wants  to  properly  represent  this  fact. 

The  curves  for  the  sets  H.’r,  Fig.  A2,  have  a  different  tip  behavior; 
they  converge  toward  a  finite  limit  curve  which  has  a  vertical  tip  tangent. 
This  limit  function  is  determined  by  the  initial  coefficients  of  the  limit 
set.  Namely,  for  r-*oo 
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A2 


C0-ri;  C2  -  -lj  C2jt-#>0  if  Xi2  (A3) 

The  limit  function  in  Fig.  A2  is  hence 

cj(y)  =  3<l-y2)V4  (r-*oo)  (A3a) 

On  the  other  hand,  the  tip  behavior  of  the  curves  in  Fig.  A1  is  produced 
by  the  tail  behavior  of  the  sets  Fjj^: 


oo 

*!»~2 

Mr 


^r/KV 


(A2a) 


We  can  use  this  fact  to  determine,  by  means  of 
of  the  curves  for  X"^  and  X'^  in  Fig.  A2.  The 


Eq,  (A2),  the  tip  behavior 
result  is 


CE(y)E=l|l-2  =  1  £ylog(2/y)  +  ... 

CE(y)EaK-3  =  'iylog(2/y)  +  ...  (A4) 


see  (1.33).  There  would  be  no  point  in  pursuing  this  comparative  procedure 
further.  Already  for  r=4  it  yields  a  leading  term,  O(y^)  ,  which  is  over¬ 
shadowed  by  the  limit  function,  Eq.  (A3). 

Also  plotted  in  Fig.  A2  are  the  curves  for  the  logarithmic  sets 
and  x'^log^Jt.  The  former  has  C2  =  -1  and  Cq  1.  The  contributions 
of  the  first  two  Legendre  polynomials  to  its  function  C^Cy)  thus  correspond 
to  those  of  the  sets  X'r,  Accordingly,  its  curve  fits  well  between  those 
for  X-3  and  On  the  other  hand,  in  the  case  of  the  set  x-^log2^we  have 

Co  =  0  and  Cq  small.  Therefore,  its  curve  remains  close  to  the  zero  line 
also  over  the  main  part  of  the  wing  span. 

By  combining  the  curves  for  X"^,  X”^  and  X’^log^K  of  Fig.  A2,  we  can 
form  the  curve  for  the  span  loading  due  to  the  general  term  of  Eq.  (44). 
However,  taken  by  itself  this  span  loading  would  not  be  very  meaningful. 

Of  the  three  conditions  that  the  terms  of  orders  s'-^log  s,  s"3  and  s"4log  s 
in  ws  must  be  zero,  it  does  fulfill  the  first  and  the  last  but  not  the 
second.  To  fulfill  the  second  condition  also,  we  add  the  set  X"^  (other 
sets  could  be  used  instead,  e.g.  the  set  Eq.  (A3)).  We  then  get  a  refined 
form  of  the  "basic  solution"  (1.48). 


4  < 


A3 


First,  wc  reformulate  the  latter  to  conform  to  the  modification 
Eq.  (47).  Thus 


It, basic*  K(Akn) 


+  V** 


<A5) 


The  second  conation  ir  fulfilled  by  setting 

a4  “  5/[  8j(2)  +  2^(3)]  =  0.3212627.. 


To  fulfill  also  the  third  condition,  we  have  to  add  the  third  general  term: 


4 

basic**  '*(4*4-1) 


31og2* 

16  TV 


(A6) 


For  this  the  second  condition  yields 


5  +  3  [4^M2)  +  ]T/K3)J/8TT2 
8f(2)  +  2}(3) 


0. 3412738. . 


whe  re 


J 


a 


(r)  = 


^lo£_K 


a:  3  specifically 

^4^(2)  -  1.989280.  ;  J*^(3)  *  0.239747.. 

The  span  loading  curves  for  these  two  basic  sets  are  shown  in  Fig.  A3. 
The  vertical  scale  of  this  figure  is  stretched.  Repeated  for  comparison  is 
the  curve  r=2  from  Fig.  A2;  this  curve  cuts  over  the  basic  curves  near  the 
tip  (compare  Fig.  4  of  Ref.  1)  because  its  log  term  in  Eq.  (A5)  is  too  large 
by  a  factor  2  for  the  first  of  the  above  three  conditions  to  be  fulfilled. 

Also  shown  in  Fig,  A3  is  the  planar  solution;  it  also  is  here  nor¬ 
malized  to  G*(l)  =  1.  This  solution  approaches  the  basic  solutions  asymptot¬ 
ically  at  the  tip,  as,  indeed,  must  any  solution  (for  any  downwash  w(y);. 


4ft 


A4 


That  tin  planar  solution  agrees  well  with  the  basic  solutions  over  the  inner 
part  of  the  wing  span  also  is  fortuitous.  In  this  latter  respect  other 
solutions  are  quite  different,  see  Fig.  ?. 


To  the  illustrations  of  elementary  span  loadings.  Figs.  A1  to  A3, 
Fig.  A4  adds  a  few  illustrations  of  corresponding  elementary  dovmwash  dis¬ 
tributions.  Denoted  by  Is  the  curve  due  to  alone;  this  curve  does 

not  fulfill  any  one  of  the  three  conditions  and  goes  to  -  oo  at  the  tip. 
Tha  curve  W£  -  W3/4,  due  to  fulfills  the  first  condition  and 


hence  approaches  the  tip  more  smoothly  than  the  first.  For  the  last  two 
curves,  a  third  term  is  added  to  fulfill  the  second  condition  also;  thus 
here  w(l)  is  finite  (but  w'(l)  is  infinite).  In  one  of  the  two  curves  the 
added  set  is  the  infinite  set  n "4  (that  is,  this  is  the  ‘Wale  solution" 
(1.43));  in  the  case  of  the  other  curve,  the  finite  set  Eq.  (A3)  is  added.  One 
sees  that  the  "basic  solution"  approaches  the  case  of  a  constant  downwash 
reasonably  well.  On  the  other  hand,  the  use  of  the  set  Eq,  (A3)  introduces  a 
curvature  corresponding  to  that  of  the  limit  function  Eq.  (A3a)  . 
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Fig,  AA.  Elementary  Downwash  Distributions 
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APPENDIX  B 

SUMMATION  FORMULAS 


In  the  analysis  of  this  paper  the  sum 


S 


E 


•± 


E 


K 


(Bl) 


of  an  infinite  set  EH  is  frequently  required.  For  example,  one  needs  Sjj  to 
determine  the  leading  coefficient  Eq  by  means  of  Eq.  (2) .  Since  the 
are  progressions,  term-by-term  summation  is  not  a  practical  proposition 
(the  rate  of  convergence  slows  down  indefinitely  as  the  summation  proceeds). 

Numerical  tables  (e.g.,  of  the  J -function  or  the  ^-function)  arc 
sometimes  available  from  which  SE  can  be  read  either  directly  or  after  some 
transformation.  This  is  not  the  case,  for  example,  if  E*  contains  logit. 
Still,  high  numerical  accuracy  is  obtained  by  means  of  the  formula,  obtained 
from  25.4.1  of  Ref.  (4) 


The  only  requirement  here  is  that  the  function  E(x)  which  interpolates  E^ 
can  be  integrated,  either  analytically,  or  numerically  after  transformation 
of  the  integral  onto  a  finite  range. 

Because  Eq.  (B2)  is  available,  there  is  no  absolute  need  for  the 
closed  form  summation  formulas  which  we  will  derive  next.  These  formulas 
are  of  considerable  mathematical  interest,  however,  because  they  appear  to 
be  the  first  of  their  kind,  namely,  the  first  to  involve  the  logarithm. 

Consider  the  sets  E  arc  si-ven  in  Eqs .  (33,  33a).  From 

Eq.  (15a)  follows 


co 

j^*c24+l,E  ‘  0  <“3) 


except  if  r=2.  Performing  this  summation  for  r  even,  one  obtains  the 
familiar  closed  form  expression  for  J(2n),  see  23.2.16  of  Ref.  (4  ).  No 
corresponding  expression  for  *J(2n+l)  seems  to  exist. 
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B2 


Performing  the  eunmation  Eq.  (B3)  for  r  odd,  we  obtain  a  new  kind 
of  result.  Making  use  of  23.2.20  of  Ref.  (4  )  and  of  some  simple  trans¬ 
formations,  we  find  for  i”3 


(JjO)  +  12JJ(2)log  2]/16 


(B4a) 


Correspondingly  for  r=5 


[3iy<5)  -  6J(2)J(3)  +  60J(4) log  2J/64 


(B4b) 


The  set  Ln  here  is  given  by  Eq.  (3la) .  It  contains  log  n  and  is  closely 
related  to  thc^-function. 

The  procedure  is  readily  continued  for  r”7,  9..  .  Note  that  the 
resulting  formulas,  Eqs.  (B4),  are  homogeneous  of  order  r  if  one  counts  as 
follows:  one  counts  the  arguments  of  the J -functions,  counts  the  power  of 
the  denominator  on  the  left,  counts  Ln  as  1,  and  notices  that  log  2  m  (1), 
see  23.2.19  of  Ref.  (4  ). 


The  formulas  Eqs.  (B4)  are  related  to  the  sets  Cjy+j.  E  through  the 
form  of  the  denominator  on  the  left.  In  order  to  obtain  formulas  which  are 


related  to  the  elements  Hfr  of  02^,  one  reverses  the  process,  setting 
■  (1 /£*■),  applies  S51  to  obtain  the  correspond  ing  set e  C2^  £,  and  makes 
use  of  Eq.  (2).  For  t  odd  this  leads  again  to  the  familiar  expression  for 
l(2n) .  For  t=4  and  t"6  one  finds 


_oo  L 

"  7J(3)/4 


X  -f  =  r31T(5)  "  1^y<2>Jr<3)J/4  (B5) 

n=l  tF  *“ 


Again  the  process  can  be  continued  ad  libitum. 

Particularly  interesting  because  of  its  simplicity  is  the  first 
formula  of  Eq.  (B5) .  This  formula  relates  directly  J(3),  the  sum  over  n“3, 
and  the  sura  over  n-zLn. 
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B3 


The  second  formula  of  Eq.  (B5)  is  used  in  this  paper  to  calculate 
the  coefficient  E0  of  its  set  Ejj,  see  Eq.  (C3).  For  this  purpose  it  is  very 
convenient.  On  the  other  hand,  our  new  formulas  may  not  be  more  convenient 
than  the  numerical  approach  Eq.  (B2)  if  the  given  set  EK  contain  log  H 
directly.  An  example  would  be  the  slowly  converging  sum 


T  22S_n  =  0.9375482. . 
n=2  n 


which  one  may  calculate  either  way. 

Note  that  the  two  sequences  of  new  summation  formulas,  Eqs.  (B4,5), 
arise  solely  from  the  series  presentation  Eq.  (9)  for  the  pressure  function 
p.  No  statement  about  the  downwash  w  is  involved,  not  even  a  statement  that 
w  should  be  physically  meaningful.  Rather,  the  Kutta  condition  Pt.e.=®  in 
Eq.  (9)  leads  directly  to  Eqs.  (2) , (13) , (14)  and  (15)  and  thence  to  Eqs. 

(B4, 5) .  The  salient  fact  is  that  the  Kutta  condition  enforces  duality 
between  the  two  sets  C2*  and  C2^+i. 
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APPENDIX  C 


ELEMENTARY  SUMS  CQ^E 


We  have 

oo 

C0,E  =  ‘  S  C2M,E 


for  the  leading  constant  of  the  elementary  set  EK  because  of  Eq.  (2). 
Thus  Cq  -  f(r)  for  the  elementary  sets  -C2*.  =  K'r,  with  If  the  Riemann 
function,  Table  2.3.  of  Ref.  (4).  For  the  remaining  sets  E^which  are 
used  to  construct  Tables  I  to  VI,  we  have 


(Cl) 


(C2) 


(C3) 
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APPENDIX  D 


PRESSURE  DISTRIBUTIONS 


The 

defined  by 


solutions  jij  Fig.  7 ,  can  be  considered  as  elementary  solutions, 
their  downwash  coefficients 


rl/(s+l)  (2s+l)  for  s=i 

l  0  "  s^i 


<D1) 


see  Eqs.  (57) and  (17).  A  given  downwash  w(y)  and  its  solution  can  be  built 
up  from  a  sufficient  number  of  such  elements. 

If  we  disregard  the  operator  (as  in  Eq.  (23))  we  have  ws  5# 
for  A-s.  Hence,  if  we  define 


1/(A+1)(2A+1) 

0 


for  A  *  A* 

"  A  *K* 


(D2) 


we  have  alternate  elements  A*  (=  0,1,2...)  which  are  somewhat  related  to  the 
elements  i.  They  also  can  be  used  to  build  up  any  desired  solution.  It  is 
of  some  interest  to  compare  the  two  types  of  elements. 

The  downwash  distributions  w^(y)  were  shown  in  Fig.  7;  Fig.  D1  shows 
the  w^*(y) .  The  major  difference  is  that  the  latter  all  turn  down  to  -  oo 
at  the  tip.  The  reason  for  this  is  clear  from  Eq.  (13):  the  -clving  set 
of  wA*  contains  only  even-numbered  rational  elementary  se'  = 

It  is  evident  from  Fig.  D1  that,  by  superposing  a  sufficient  r.urr.nar®  of 
elements  A*,  one  can  approximate  a  given  w(y)  closely  enough  within  any 
given  inner  range/y|<y0<l.  However,  in  order  to  achieve  w(l)  finite,  one 
has  to  represent  the' set  and  this  requires  an  infinite  number  of 

elements  X*. 

The  purpose  of  the  present  appendix  is  to  discuss  pressure  functions 
p.  We  are  interested  in  cases  where  w(y)  is  wavy,  and  are  interested  in 
the  nature  of  the  pressure  distribution  over  the  inner  part  of  the  wing  span. 


4)  In  some  respect,  this  corresponds  to  the  procedure  of  collocation  analyses. 


From  this  point  of  view,  the  elements  i  and  \*  are  equivalent.  Figures 
D2  and  D3  show  the  cases  X*  =  1  and  A*  =  3  (instead  of  the  corresponding 
cases  i  c  1,3  simply  because  the  present  figures  had  already  been  prepared). 

Like  Fig.  2,  Figs.  D2  and  D3  ore  relief  diagrams,  however,  the  right 
half-wing  is  stretched  into  a  rectangle;  furthermore,  the  rectangle  is  cut 
at  y  =  0.4  for  bcLter  visibility  of  p.  (The  wing  tip  point  is  stretched 
into  a  straight  line,  a  "chord"  y  =  1.0;  thus  the  linear  distribution  of  p 
over  the  front  part  of  this  "chord"  is  shown.)  A  vertical  plane  through 
the  leading  edge  of  the  rectangle  is  drawn;  it  shows  horizontal  lines  of 
constant  p  (the  scale  of  p  is  arbitrary)  and  vertical  lines  of  constant  y. 

The  left  half-wing  is  not  shown  in  Fig.  D2.  It  is  shown,  unstretched, 
in  Fig.  D3,  and  over  it  the  wavy  incidence  distribution  w(y)  is  indicated. 

The  pressure  distributions  Figs.  D2, 3  lead  to  a  somewhat  unexpe.cted 
observation  of  technical  interest:  while  the  waviness  of  the  incidence 
distribution  is  well  reflected  in  p  along  the  leading  edge,  this  waviness 
is  rapidly  damped  out  along  the  chord  and  already  at  midchord  is  no  longer 
recognizable.  Over  the  rear  part  of  the  wing,  the  pressure  distribution 
p  no  longer  reflects  the  waviness  of  w(y) .  This  stabilizing  effect  (shown 
here  on  an  idealized  model)  must  be  generally  significant  for  the  influence 
of  disturbances  of  short  characteristic  spanwise  extension.  It  is,  of 
course,  also  the  mechanism  which  makes  the  wave  amplitudes  in  ^  of  Fig.  7 
smaller  than  the  corresponding  amplitudes  in  wi. 

Another  way  of  describing  the  mechanism  which  leads  to  Fig.  D3  in 
particular  is  to  say  that  already  the  vortices  which  are  shed  by  the  front 
wing  create  a  downwash  field  which  is  almost  identical  with  that  of  the  back 
part  of  the  wing,  such  that  at  the  back  part  almost  no  new  vorticity  is 
created.  This  mechanism,  respectively  the  chordwise  distributions  of  the 
lift  which  are  shown  ir  Fig.  D3,  are  distinctly  different  from  the  assumptions 
of  lifting  line  theory.  Accordingly,  lifting  line  theory  cannot  be  expected 
co  perform  adequately  for  the  downwash  wj.  In  Fig.  D4,  the  curve  3  of 
Fig.  7  is  repeated  (one  half -wing  only,  larger  vertical  scale),  and  the 
result  of  lifting  line  theory  is  also  shown.  Comparison  of  the  two  results 
reminds  one  of  the  fact  that  lifting  line  predicts  the  lift  coefficient  of 
the  planar  slender  elliptic  wing  too  large  by  the  factor  2.  According  to 
Fig.  D4,  this  amplification  factor  2  applies  also  to  wings  which  have  a 
wavy  distribution  of  the  wing  incidence,  supposing  of  course  that  the  wave 
lengths  are  sufficiently  small. 
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Fig.  D4.  Comparison  of  Lifting  Line  Result  with  Exact  Solution 


